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Abstract 

We develop a Hamiltonian theory for 2D soliton equations. In particular, we identify 
the spaces of doubly periodic operators on which a full hierarchy of commuting flows can 
be introduced, and show that these flows are Hamiltonian with respect to a universal 
symplectic form u = |ReSoo < \E'o5LA5\E'o > d/c. We also construct other higher order 
symplectic forms and compare our formalism with the case of ID solitons. Restricted 
to spaces of finite-gap solitons, the universal symplectic form agrees with the symplectic 
forms which have recently appeared in non-linear WKB theory, topological field theory, 
and Seiberg-Witten theories. We take the opportunity to survey some developments in 
these areas where symplectic forms have played a major role. 
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I. INTRODUCTION 



There is increasing evidence that symplectic structures for solitons may provide a 
unifying thread to many seemingly unrelated developments in geometry and physics. In 
soliton theory, the space of finite-gap solutions to the equation [dy — dt — A] = is a 
space A4g{n,m) of punctured Riemann surfaces F and pair of Abelian integrals E and 
Q with poles of order less than n and m at the punctures. The fibration over A4g{n, m) 
with F as fiber carries a natural meromorphic one-form, namely dA = QdE. It is a 
remarkable and still mysterious fact that the form dA is actually central to several theories 
with very distinct goals and origins. These include the non-linear WKB (or Whitham) 
theory [22] [23] [26] [36] [38] , two-dimensional topological models [11] [12] [13], and Seiberg- 
Witten exact solutions of N=2 supersymmetric gauge theories [18] [40] [52] [53]. The form 
dA can be viewed as a precursor of a symplectic structure. Indeed, it can be extended as a 
1-form XliLi d^i^i) on the fibration over Aig{n, m) with fiber a symmetric ^rth-power of F. 
Its differential uj becomes single- valued when restricted to a suitable (/-dimensional leaf of 
a canonical foliation on A^p(n, m), and defines a symplectic form [39]. Earlier special cases 
of this type of construction were pioneered by Novikov and Veselov [50] in the context of 
hyperelliptic surfaces and ID solitons, and by Seiberg, Witten, and Donagi [18] [53] in the 
context of N=2 SUSY gauge theories. 

The goal of this paper is twofold. Our first and primary objective is to construct the 
foundations of a Hamiltonian theory of 2D solitons. 

• For this, we provide an improved formulation of 2D hierarchies, since the classical 
formulations (e.g. Sato [54]) are less pliable than in the ID case, and inadequate for 
our purposes. In particular, the new formulation allows us to identify suitable spaces 
C{h) of doubly periodic operators on which a full hierarchy of commuting flows dmL = 
dyAm -\- [Ami L] can be introduced; 

• We can then define a universal symplectic form uj on these spaces C{h) by 

a; = ^Resoo < *o^-^ /\ > dA; (1.1) 

where and '^q are the formal Bloch and dual Bloch functions for L. This form had been 
shown in [39] to restrict to the geometric symplectic form (^X]i=i dA(2;i) when finite-gap 
solitons are imbedded in the space of doubly periodic operators. Here we show that it is a 
symplectic form in its own right on £(6), and that with respect to this form, the hierarchy 
of 2D flows is Hamiltonian. Their Hamiltonians are shown to be nHn+m, where Hg are 
the coefficients of the expansion of the quasi-momentum in terms of the quasi-energy. 

• Our formalism is powerful enough to encompass many diverse symplectic structures 
for ID solitons. For example, uj reduces to the Gardner-Faddeev-Zakharov symplectic 
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structure for KdV, while its natural modifications for y-independent equations (see (2.71) 
and (2.73) below), reproduce the infinite set of Gelfand-Dickey as well as Adler-Magri 
symplectic structures. 

• The symplectic form (1.1) is algebraic in nature. However, it suggests new higher 
symplectic forms. 



which are well-defined only on certain spaces of operators with suitable growth or ergod- 
icity conditions. For Lax equations dt^L = [Am,L], these higher symplectic forms have a 
remarkable interpretation: they are forms with respect to which the eigenvalues of A^^, 
suitably averaged, can serve as Hamiltonians, just as the eigenvalues of L are Hamiltoni- 
ans with respect to the basic symplectic structure (1.1). It would be very interesting to 
understand these new forms in an analytic theory of solitons. 

Our second objective is to take this opportunity to provide a unified survey of some 
developments where the form dA (or its associated symplectic form uj) played a central 
role. Thus dA emerges as the generating function for the Whitham hierarchy, and its 
coefficients and periods are Whitham times (Chapter IV). The same coefficients are de- 
formation parameters of topological Landau-Ginzburg models in two dimensions (Chapter 
V), while for N=2 SUSY four-dimensional gauge theories, the periods of dA generate the 
lattice of Bogomolny-Prasad-Sommerfeld states (Chapter VI). Together with dA, another 
notion, that of a prepotential JF, emerges repeatedly, albeit under different guises. In non- 
linear WKB methods, JF is the exponential of the r-function of the Whitham hierarchy. In 
topological Landau-Ginzburg models, it is the free energy. In N=2 supersymmetric gauge 
theories, it is the prepotential of the Wilson effective action. It is an unsolved, but clearly 
very important problem, to determine whether these coincidences can be explained from 
first principles. 



Solitons arose originally in the study of shallow water waves. Since then, the notion 
of soliton equations has widened considerably. It embraces now a wide class of non-linear 
partial differential equations, which all share the characteristic feature of being expressible 
as a compatibility condition for an auxiliary pair of linear differential equations. This is 
the viewpoint we also adopt in this paper. Thus the equations of interest to us are of the 
form 



mo — 




(1.2) 



II. HAMILTONIAN THEORY OF 2D SOLITON EQUATIONS 



[dy-L,dt-A]=0, 



(2.1) 
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where the unknown functions {ui{x,y,t)}f^Q, {vj{x,y,t)}jLQ are the N x N matrix coef- 
ficients of the ordinary differential operators 

n m 

L = Y,Ui{x,y,t)dl A = Y,'t^^{^.y.t)di (2.2) 

A prehminary classification of equations of the form (2.1) is by the orders n, m of the 
operators L and A, and by the dimension N of the square matrices Ui{x,y,t), Vj{x,y,t). 
In what follows, we assume that the leading coefficients of L and A are constant diagonal 
matrices u^^ = u^5af3, = v^Sap- Under this assumption, the equation (2.1) is invari- 
ant under the gauge transformations L,^ — > L' = g{x)~^Lg{x), A' = g{x)Ag{x)~^ where 
g{x) is a diagonal matrix. We fix the gauge by the condition ■u^'^ = d^^Ua^ u'^-i = 0. 
We shall refer to (2.1) as a zero curvature or 2D soliton equation. The ID soliton equa- 
tion corresponds to the special case of y-independent operators L and A. In this case the 
equation (2.1) reduces to a Lax equation = [A,L]. 

A. Difficulties in a Hamiltonian Theory of 2D Solitons 

The Hamiltonian theory of ID solitons is a rich subject which has been developed 
extensively over the years [10] [25]. However, much less is known about the 2D case. 
We illustrate the differences between ID and 2D equations in the basic example of the 
hierarchies for the Korteweg-deVries (KdV) 

Ut - ^ud^u + ^dlu = (2.3) 

and the Kadomtsev-Petviashvili (KP) equations 

3 3 3 

-Uyy = d^{ut - -ud^u + -d^u). (2.4) 

The KP equation arises from the choice N = 1, n = 2, m = 3, and L = 9^ + u, A = 
dl + ^udx + vo in (2.1). We obtain in this way the system 

dxvo = ^dlu + ^Uy, vo,y = ^d^u + ^d^Uy - ^udxu (2.5) 

which is equivalent to (2.4) (up to an (a;, y)-independent additive term in vq, which does 
not affect the commutator [dy — L,dt — A]). Taking L and A independent of y gives the 
KdV equation. 

The basic mechanism behind this construction is that the zero curvature equation 
actually determines A in terms of L. This remains the case for the ID Lax equation 
Lt = [A, L] even when A is taken to be of arbitrarily high order m, but not for the 2D zero 
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curvature equation L-t — Ay = [A, L] . The point is that [A, L] is a differential operator of 
order m + 1. The Lax equation requires that it be in fact of order 0, while the 2D zero 
curvature equation requires only that it be of order < m — 1. The order constraint is 
quite powerful. Expressed as differential constraints on the coefficients of A, it implies 
readily that the space of such A's for fixed L is of dimension m. An explicit basis can 
be obtained by the Gelfand-Dickey construction [10] [27], which we present for a general 
operator L of order n. Let a pseudo-differential operator of order n be a formal Laurent 
series Yl^=-oo '^i^x with and satisfying the identities 

oo 

d^U = Uda, + U', d-^U = J2i-y'^^^d-'-\ 

Then there exists a unique pseudo-differential operator L^/" of order 1 satisfying (L^/"^)" = 
L. Evidently, the coefficients of L^/" are differential polynomials in the coefficients of L. 
For example, for L = d^ + w, we find L^/^ = dx + \ud~^ — \u'd~^ + • • • . We set 

where the first term on the right hand side is the differential part of the pseudo-differential 
operator L*/"^, and the second term on the right hand side is of order < —1. Then 
[L, lX""] = [L, U/""] - [L, L!_/"] = -[L, L!/'']. Since the commutator [L, L!/"] is of order at 
most n — 2, this shows that the differential operators L^Jj^ provide the desired basis. 

Associated to L are then an infinite hierarchy of fiows, obtained by introducing "times" 
ill' ■ ■ 1 tm, ■ ■ •) and considering the evolutions of L = Yl7=o ^(^' ^i; ' ' ' > ^n)^^ defined by 

dmL=[L,L'^/'']. (2.6) 

where we have denoted by dm the partial derivative with respect to the time tm- A key 
property of these fiows is their commutativity, i.e. 

[d. _ L!/-, dj - L'l''] = 0. (2.7) 

To see this, we note first that if L evolves according to a fiow dtL = [L, A\, then evolves 
according to a^L" = [L", A]. Thus we have diU/"" - [L^/", L^/"], a^L^/" = [L^^", L*/'^], 
and the left hand side of (2.7) can be rewritten as -[lX"" , U/''] + [L^'', D/''] + [L%'' , 
If we replace VJJ^ by L*/" — Vly' and U^"^ — L^J^, all terms cancel, except for 

[L!^", lI/"']. This term is however pseudo-differential, of order < —3, and cannot occur in 
the left hand side of (2.7), which is manifestly a differential operator. 
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The flows (2.6) are known to be Hamiltonian with respect to an infinite number of 
symplectic structures with different Hamiltonians. For example, the KdV equation itseff 
can be rewritten in two Hamiltonian forms 

SH SH' 

ut = dx — , ut = {dl + 2{udj; + da;u))-^ 

where the skew-symmetric operators K = d^, K' = + 2{udx + d^u) correspond to two 
different symplectic structures (called respectively the Gardner- Faddeev-Zakharov [10] [25] 
and Adler-Magri structures [1][41]), and H — — ^u^, H' — are the corresponding 
Hamiltonians. 

The situation for the 2D zero curvature equation is much less simple, since the argu- 
ments narrowing A to an m-dimensional space of operators break down. Although formally, 
we may still introduce the KP hierarchy as dmL — dyAm + [Am, L], with Am an operator 
of order m which should also be viewed as an unknown, this is not a closed system of 
equations for the coefficients of L, as it was in the case of the Lax equation. Another way, 
due to Sato [54] , is to introduce the KP hierarchy as a system of commuting ffows 

dml^=[L^M (2.8) 

on the coefficients {vi{x,t))lZ'^ of a pseudo-differential operator L 

CO 

l. = d + ^Vi{x,t)d-\ 

i=l 

In this form, the KP hierarchy can be viewed as a completely integrable Hamiltonian 
system. However, it now involves an infinite number of functions Vi, and its relation to the 
original KP equation (which is an equation for a single function of two variables {x,y)) 
requires additional assumptions. 

A. Quasi-Energy and Quasi-Momentum 

Our first main task is then to identify the space of differential operators with periodic 
coefficients on which the KP equation and its higher order analogues can be considered as 
completely integrable Hamiltonian systems. Our approach actually applies systematically 
to general 2D soliton equations. We present these results at the end of this section, and 
concentrate for the moment on the simplest case of a differential operator L of order n. We 
begin with the construction of the formal Bloch eigenfunction for two-dimensional linear 
operators with periodic coefficients. 

Theorem 1. Let L be an arbitrary linear differential operator of order n with doubly 
periodic coefficients 

n-2 

Ui{x + l,y) = Ui{x,y+1) = Ui{x,y) (2.9) 



6 



Then there exists a unique formal solution ^o{x,y; k) of the equation 



{dy-L)^o{x,y;k) = 0, (2.10) 

which satisfies the following properties 

(i) ^o{x, y; k) has the form ^o{x, y; /c) = (1 + U^, y)k-') e ^+^-0' ; 

(ii) ^'o(0,0;fc) = l; 

(iii) '^q{x, y; k) is a Block function with respect to the variable x, i.e., 

'^o{x+l,y,k) = wi{k)'^o{x,y,k), wi{k) = e'' (2.11) 

The formal solution ^o{x, y, k) is then also a Bloch function with respect to the variable 
y with a Bloch multiplier W2{k) 

mo{x, y + 1; /c) = W2{k)^o{x, y; k), W2{k) = (1 + ^ j,A;-«)e('="+Sr="o (2.12) 

s=l 



Proof. To simplify the notation, we begin with the proof in the case of n = 2, with 
L — d^-\- u{x, y). The formal solution has then the form 

Substituting this formal expansion in the equation [dy — L)^o{x,y; k) = gives the fol- 
lowing equations for the coefficients 

dyis + dyBois = K+i + C + <s. (2.13) 

(Here and henceforth, we also denote derivatives in x by primes.) These equations are 
solved recursively by the formula: 

6+1 {x,y) = cs+i {y) + ^s+i (x, y) , 

^s+ii^^y) = \j^ {dyUx:y)-C{X:y) + {dyBo - u{x,y))Ux:y))<^5: (2.14) 

where Cs{y) are arbitrary functions of the variable y with the only requirement that Cs(0) = 
0, which is dictated by (ii). 

Our next step is to show by induction that the Bloch property (i), which is equivalent 
to the periodicity condition 

es(a;+l,y) = e.(a;,y), (2.15) 
7 



uniquely defines the functions Cs{y). Assume then that is known and satisfies 

the condition that the corresponding function ^°(a;,y) is periodic. The first step of the 
induction, namely the periodicity in x of ^i(a;, y), requires that dyBo{y) = +u, in view 
of (2.13) and the fact that = 1- Thus we fix Bo{y) to be 

Bo{y)= r [ u{x',y')dx'dy' (2.16) 
Jo Jo 

The choice of the function Cg{y) does not afi^ect the periodicity property of ^s{x,y), but 
it does affect the periodicity in x of the function ^s+ii^^v)- order to make ^s+ii^^v) 
periodic, the function Cs{y) should satisfy the equation 

dyCs = - f {dyes{x,y) - (e°)"(^,y) + {dyBo - u{x,y))es{x,y))dx (2.17) 

where < • >x denotes averaging in x over the interval [0,1]. Together with the initial 
condition Cs(0) = 0, this defines Cg uniquely 

cs = - fdy' ['{dy^A^.y') - {0"{^:y') + {dyBo - u{x,y'))CAx.y'))dx, (2.18) 
Jo Jo 

and completes the induction step. 

We can now establish (iii), also by induction. Assume that ^^-i satisfies the relation 

s-1 

^s-i{x,y+l) -^s_i{x,y) = ^Ji^s-i-i{x,y). (2.19) 

i=l 

where Ji, . . . , Js_i are constants. Then (2.14) implies that 

s — 1 s 

Cs{x, y + 1)- Cs{x, y)^^ Ma-i{x, y) + Js^^ Ms-i{x, y), (2.20) 

i=i i=i 

with 

s-1 

Js = Cs{y + 1) - Cs{y) - ^ Cs-i{y). 
We claim that Jg is actually constant. In fact, it follows from (2.14) and (2.19) that 

s-1 s-1 

cAx, y + 1)- ^s (^, y)-J2 y)-J2 Mis-^{x, y) - c^-M). 

i=i i=i 
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Thus (2.18) implies that 

dycsiy + 1) - dycsiy) = - £ f\dyes-i - iO" + {dyBo - «)e°)da; 

s-l s-1 



i=\ 1=1 



In particular the derivative of Jg vanishes. This proves Theorem 1 when n — 2. 

The proof can be easily adapted to the case of general n. Let (x, y) be the coefficients 
of the formal series 

= w^^) *^)*-')«'^- 

Then [dy — L)^'o = is equivalent to the system of equations 

n— 1 n— 2 i n—2 

E + E E (^r^e°+0 = dy^'s + E ^^■^?+- s>-n + 2, (2.21) 

i=o i=0 i=o j=-s 

where we have assumed that = for s < 0, and set hj — bj{y) to be the coefficients of 
the series 

^yM^^k-+ y bMk^. (2.22) 

These equations are of the form 

nd^^s+n-i = b-s + FsiC, fcs"), (2.23) 

Here F is a linear combination of the s' < s + n — 2 and their derivatives, with 
coefficients which are themselves linear in Ui and bg", s" < s (c.f. (2.13)). The equations 
define recursively b-s{y) and The coefficient b-s follows from the periodicity in x 

of (2.23) 

6_,(y) = - / Fs^CMAx, (2.24) 
and the coefficient Cs+n-i follows in turn 

nC+n-i = b-s{y)x + f FsiCMdx'. (2.25) 
We can now integrate (2.22) and find that has the form (i) with 

Bi{y)^ r bi{y), i = 0,...,n-2. (2.26) 
9 



We obtain at the same time the Bloch property for with respect to y with the Bloch 
multipher 

»1 n-2 

W2{k)= {k''+ T h{y)k')dy. (2.27) 
The proof of Theorem 1 is complete. 

B. Basic Constraints 

In Theorem 1, we have chosen the simple form wi (k) = for the Bloch multiplier in x. 
If we view a; as a "space" variable, this identifies the spectral parameter k with the "quasi- 
momentum" (up to a factor of z = -\/— T) . The variable y can then be interpreted as a time 
variable in the Schrodinger-like equation [dy — L)'^q = 0. This identifies (again up to a 
factor of i) the logarithm of the second Bloch multiplier W2{k) with the quasi-energy E{k). 
Alternatively, we may change spectral parameters, and introduce the spectral parameter 
K as well as the coefficients (£^i)^_^^.2 of the expansion in k of the quasi-energy by 

oo 

K"" = E{k) = \ogW2{k) = k'' + Sik-' 

i=-n+2 

K = k + 0{k-^) (2.28) 

We observe that a change of spectral parameter of the form k ^ K = k+0{k~^) transforms 
a formal series in k of the form (i) in Theorem 1 into another formal series of the same 
form in K. With K as spectral parameter, the second Bloch multiplier W2{K) reduces 
to W2{K) — exp(iir"), but the first Bloch multiplier wi{K) becomes non-trivial. The 
coefficients {Hs)'^i of the expansion in K of its logarithm 

oo 

wi{K) = e''^^\ k{K) = K + J2HsK-'. (2.29) 

will play an important role in the sequel. 

The coefficients Si and Hg of the expansions (2.28) and (2.29) of the spectral param- 
eters are uniquely defined by the coefficients {ui)^~Q of the operator L, and hence can be 
considered as functionals on the space of periodic operators. We now restrict ourselves to 
the subspace C{b) of operators L satisfying the constraints 

h{y) = 6^, < z < n - 2, (2.30) 

where bi{y) = dBi{y)/dy, with the functions Bi{y) defined by (i) in Theorem 1, and 
b = {bo, • • • , bn-2) are (n — 1) fixed constants. On the space C{b), the essential singularity 
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in the second Bloch multiplier W2{k) simplifies to k'^ + Yl7=o (^•^- (2-12)). Comparing 
with (2.28), we see that the constraints (2.30) fix the values of the first (n — 1) functionals 
£-i = hi- This is in turn clearly equivalent to fixing the values of the first {n — 1) functionals 
Hi, . . Hn-i in (2.29). We claim that these constraints can also be expressed under the 
form ^ 

/ hi{uj)dx = Hg = const, z = 1, . . . , n — 1, (2.31) 

Jo 

where hs(uj), j > (n — i — 1), are universal difi^erential polynomials depending only on 
n. In fact, the constraints (2.30) imply that the essential singularity in the Bloch function 
^o{x, y; k) is of the form 

n-2 

kx + {k'' + bik')y ^Kx + K'^y + 0{K-^) 

i=0 

Since the expression exp (0{K~^)) contains no essential singularity and can be expanded 
as a formal series in K~^, we have shown that the formal Bloch solution \E'o(a;, y; fc) can 
be rewritten in terms of the spectral parameter K under the form 

*o(x, V^K)=(^ + Y. V)K-^ e^-+^"^ (2.32) 

Substituting (2.32) in [dy — L)'^ q{x , y; K) = gives for the first (n — 1) coefficients 
(i, . . . , Cn-i a system of ordinary differential equations 

n—l n—2 i 

E {dr'cUl) + E E (^r^co = o, . = -n + 2 . . . , o, (2.33) 

1=0 i=0 1=0 

which just coincides with the first equations defining formal eigenf unctions for ordinary 
differential operators (see [kr]). Let Ci{x,y;xo), be the solution of (2.33) with the nor- 
malization ^0 = 1 and Q{xo,y]XQ) = 0. Then the equations (2.33) define recursively 
differential polynomials hi{uj{x,y)) such that 

dxCi{x,y;xo)\x=xo = hi{uj{x,y))\x=xo- (2-34) 

The above left hand side is equal to the first coefficients of the logarithmic derivative of 
^0 ai x = xo, i.e. 

^ + E (^-C.(^, y, xo)\xo=.) K-^ = ^-^"(^-^-^^ + 0{K-). (2.35) 

^o[,x,y,J\) 

Integrating gives the first (constant) coefficients of (2.29) and establishes our claim. 
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Henceforth we will always assume that L is in the space C{h) of operators with periodic 
coefficients {ui)'^~Q satisfying either one of the equivalent constraints (2.30) or (2.31). 

C. Commuting Flows 

On the space C{h) we can now define an infinite set of mutually commuting fiows as 
follows. First, we observe that for any formal series of the form (i) in Theorem 1, there 
exists for each integer m > 1 a unique differential operator A^^ of the form 

m— 1 
3=0 

which satisfies the condition 

(Am - K'^ik))Mx, y; k) = Oik-')^^, y; k), (2.37) 

where K{k) is defined by (2.28). Indeed, this condition is equivalent to the following finite 
system of equations for the coefficients (tti,m) of A^ 

m i 

^^^^.^^^^(^^(^r'Cs+O = Cs+m, s = -m + 2,...,0. 

This system is triangular, and identifies uniquely the coefficients Ui^rn ^ differential poly- 
nomials in the first m — 1 coefficients of ^o{x, y; k). For example, we find 

Um-2,m = -mC[, Um-3,m = -mC'2 - C^Cl 

Let L be now an operator in Cih), '^^{x.y.k) its Bloch function, and A^ be the 
operators obtained by the preceding construction. Then 

[dy - L,Am]'l>o{x,y,k) = {dy - L)0{k-^)^o{x,y,k) = 0(/c"-2)e^^+^"^ (2.38) 

which implies that the operator [dy — L, A^] has order less or equal to n — 2. We set 

n-2 

[dy - L, Am] = dyAm - [L, Am] = ^ F,,m{x, y)di (2.39) 

The functions Fi^^ix, y) are uniquely defined by Uj{x, y) and can by expressed in terms of 
multiple integrals of differential polynomials in Uj. Thus they can also be expressed as 

Fi^rn{x,y) = Fi,m{u{x,y)), U= {uo,...,Un-2)- (2.40) 

12 



(2.36) 



where is a nonlocal hut exact functional of Uj{x,y). 

Theorem 2. The system of equations 

dmUi = Fi^rn{u) < — > = dyA^ + [A^, L] (2.41) 

defines an infinite set of commuting flows on the subspace C{h) of doubly periodic operators. 

Since Fj^TO and A^, are well-defined functions of L, we need only check the commu- 
tativity of the flows. For this, we need the following version of the uniqueness of Bloch 
solutions, namely, that if is a Bloch solution [^{x+ l,y;k) = wi{k)'if{x,y;k)) 

of the equation {dy — L)^ = having the form 

M^,y;k) = ( Ux,y)kA e''^+^'"+^:::'^'''^y, (2.42) 

\s=-N J 

then ^{x,y; k) must be given by 

^{x,y;k) = a{k)^o{x,y;k), (2.43) 

where a{k) is a constant Laurent series a{k) = YI'^n ct'sk~^- This is because the leading 
coefficient ^jv in any formal solution is a constant. As a consequence, if '^{x,y;k) is a 
Bloch solution, then the ratio ^{x, y; k)'^~^{0, 0; k) is also a solution. This expression has 
all the properties of, and can be identified with ^^o- Since ^^(0,0; A;) is a Laurent series, 
our assertion follows. 

Returning to the proof of Theorem 2, we observe that if ^'o is the Bloch solution to 
(2.22) then {dm — ^m)*o is also a Bloch solution to the same equation. The preceding 
uniqueness property implies 

(a^-A^)^'o = -fim(A:)*o, n^{k)^K^ + 0{k~^). (2.44) 

for a suitable {x,y) independent Laurent series Qmik)- In particular, [dm — Am,dn — 
An]'^o = 0{k~^)'^o- The last equality implies that [dm — Am,dn — An] is an ordinary 
differential operator in x of order less than zero. Therefore, it must vanish identically 

[dm-Am,dn-An]^0 (2.45) 

This establishes the commutativity of the ffows. 

Conversely, let u{x, y,t,k), t — {ti, t2, . . .) be a solution to the hierarchy (2.41). Then 
there exists a unique formal Bloch solution ^'(a;, t; K) of the equations 

{dy-L)^ = 0, {dm-Am)^ = (2.46) 
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having the form 

^{x, y, t,K)=(^l + f2 (six, y, t)K-^^ e^^+^"^+S. (2.47) 

The above expression identifies the original variables x and y with the first and the 
n-th times of the hierarchy, 

x = tu y = tn, (2.48) 

respectively. More generally, the preceding results show that for periodic operatos L, an 
equation of the zero curvature form [dy — L^dt — A\ = must be equivalent to a pencil of 
equations for the coefiicients of L only. In other words, there must exist constants Ci such 
that 

m 

A = Y, a^i (2-49) 

i=l 

and the flow is along the basic times ti = Cit, of the hierarchy. 

Finally, we point out that for all n, the equations of the corresponding hierarchy for 
L = Ln have the same form [dn — Ln, dm — Lm\ = 0, and can be considered as reductions 
of this system. Our approach to these reductions is to select two particular times which 
we treat as spatial variables, and to impose periodicity conditions in these variables. 

D. Dual Formal Bloch Solutions 

A key ingredient in our construction of symplectic structures on spaces of periodic 
operators L is the notion of dual Bloch functions ^q{x, y; k). In the one-dimensional case, 
dual Bloch functions were introduced in [7]. In our set-up, its main properties are as 
follows: 

• Let \E'o(2;, y, k) be a formal series of the form (i) in Theorem 1. Then there exists a 
unique formal series ^o(^; V'l ^) form 

**(x,2/;fc) = e-^--(^"+5:r:;''^^^)^ (^i+pcAx,y)k-'^ , (2.50) 

such that for all non-negative integers m the equalities 

Resoo{%{x,y;k)d^^o{x,y;k)dk) = 0, m = 0,l,... (2.51) 

are fulfilled; 

• If \l/o is a Bloch function with Bloch multipliers Wi{k), then \E'o(x,|/; k) is a Bloch 
function as well with inverse Bloch multipliers 

%{x + 1, y; k) = w^\k)%{x, y; k), %{x, y+l;k) = w^\k)%{x, y; k), (2.52) 
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• If ^o{x, y; k) is a solution to the equation {dy — L)'^q = 0, then the series ^* {x, y; k) 
is a solution to the adjoint equation 

%{dy-L) = Q, (2.53) 

where the action on the left of a differential operator is defined as a formal adjoint action, 
i.e. for any function /* 

{rdi) = {-d,yr. (2.54) 

To see this, we begin by noting that, although each of the factors in (2.51) has an 
essential singularity, their product is a meromorphic differential and the residue is well- 
defined. It has the form 

ReSoo (^S^r^od/C) = Cm + C + 9m{^U U-i; Cl,---, Cm-l), 

where gm is linear in in and their derivatives, s < m. The condition (2.51) defines 
then recursively as differential polynomials in ^g, s = 1, . . . , m. For example, we have 

This shows the existence and uniqueness of \l/o. Since the second statement is a direct 
corollary of the uniqueness of \E'o, we turn to the proof of the last statement. First, we 
show that if ^* {x, y; k) is a formal series 

^\x,y-k) = e-'=^-('="+E:r;^^^^)^ [ Cs{x,y)kA , (2.55) 

\s=-N } 

satisfying the equations (2.51), then there exists a unique degree N ordinary linear differ- 
ential operator D such that 

**(x,y; k) = *;(x,|/; k)D. 

Since ^^^'o satisfies the equations (2.51), we can find D satisfying the condition 

**(x, y; k) - *S(x, y; k)D = 0{k-')%{x, y; k). 

The above right hand side has the form (2.55) with < and satisfies (2.51). Evaluating 
the leading term, we find that it must vanish identically. 
Let *o be a solution of {dy - L)^'o = 0. Then 

ReSoo (dy^d^'^odk) = dyReSoo {^d^'^odk) - ReSoo (*o^^^*od/c) = 
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In particular, there exists a differential operator L such that 

Let f{x) be an arbitrary periodic function on one variable. We have 

dy < f{x)%^o >x=< f{x)%{L + L)^o) >x, 

where we have denoted as usual the average value in x of any periodic function g{x) by 
< g >x- The above left hand side is of order —1 in k. On the other hand, if L + L is not 
equal to zero and gi{x^y), < i < n — 2 are its leading coefficients, then the right hand 
side is of the form < f{x)gi{x,y) >j: /c* + This implies that < fgi >x= 0. Since 

/ was arbitrary, we conclude that ^r^ = 0, establishing the last desired property of dual 
Bloch functions. 

We conclude our discussion of dual Bloch functions with several useful remarks. The 
first is that the identity 

<%{x,y;k)^o{x,y;k)>x=l. (2.56) 

holds for any formal series ^o{x,y;k) of the form (i) in Theorem 1 and its dual Bloch 
series '^Q{x,y;k). Indeed, just as in Section II. C, we can show the existence of a unique 
pseudo-differential operator $ = 1 + Yl'^o '^ii^i y)9x^ ^^at 

^o{x, y; k) = $e'="+('="+Sr=o '^^'^'^y. (2.57) 

As in [10], this implies 

%{x,y;k) = (^e-'=^-('="-Sr=o''^'=')f) (2.58) 

More precisely, let Q = YIi'^n 1-i-i^)^x^ be a pseudo-differential operator. Then we may 
define its residue resgQ by 

resgQ = ReSoo(e~''^((5e''"'))dA; = -qi. 

The point is that, while the ring of pseudo-differential operators is not commutative, the 
residue is, after averaging 

< resd{QiQ2) >x^< resd{Q2Qi) >x ■ 

This shows that the series defined by the right hand side of (2.58) satisfies (2.51), and 
hence must coincide with ^^q. The desired identity (2.56) is now a direct consequence of 
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the two preceding identities and of the associativity of the left and right actions under 
averaging. 

Secondly, we would like to stress that, although '^q is a Bloch solution of the ajoint 
equation '^o{dy — L) = 0, its normalization is different from that used for ^^o- This 
symmetry may be restored if we introduce 

The inverse relation is then 



%{x,y; k) 



Finally, the definition of the action on the left of a differential operator adopted earlier 
implies that for any degree N differential operator 



N 



i=0 



there exist degree {N — i) differential operators D^*) such that for any pair of functions /+ 
and g the equality 

n 

holds. The set of operators D^*^ was introduced in [36]. Of particular interest is of course 
= D, and the "first descendant" of D, namely 

n 

D^^^ ^J2^w,{x)dl-\ (2.60) 

i=0 



E. The Basic Symplectic Structure 

We are now in position to introduce a symplectic structure on the space C{b) of 
periodic operators L subject to the constraints (2.30), and to show that the infinite set of 
commuting flows constructed in Theorem 2 are Hamiltonian. 

The main ingredients are the one- forms 6L and S'^q. The one-form 6L is given by 

n-2 

5L = ^Su.dl., 

i=0 
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and can be viewed as an operator-valued one-form on the space of operators L = d'^ + 
Yl7=o '^i^x- Similarly, the coefficients of the series are explicit integro-differential poly- 
nomials in Ui- Thus can be viewed as a one- form on the space of operators with values 
in the space of formal series. More concretely, we can write 

(oo \ / oo \ 

The coefficients 5^s (or <^Cs ) can be found from the variations of the formulae (2.24), (2.25) 
for or recursively from the equation 

{dy- 1)5-^0 = {5^0- (2.61) 
Let /(a;, y) be a function of the variables x and y. We denote its mean value by 

<f>= [ [ f{x,y)dxdy. 
Jo Jo 



Theorem 3. (a) The formula 

uj=^ ReSoo < %SL A 5*0 > dk, (2.62) 

defines a symplectic form, i.e., a closed non-degenerate two-form on the space C{h) of 
operators L with doubly periodic coefficients, (b) The form co is actually independent of 
the normalization point {xq — 0,yo = 0) for the formal Bloch solution ^o{x, y; k). (c) The 
flows (2.41) are Hamiltonian with respect to this form, with the Hamiltonians niy^+„(w) 
defined by (2.29). 

Proof. We require the following formula, which is a generalization of the well-known 
expression for the variation of energy for one-dimensional operators. Let E{k) be the 
quasi-encrgy which is defined by (2.28). Its coefficients are nonlocal functionals on the 
space C{b) of periodic functions Ui{x,y) subject to the constraints (2.30). Then we have 

dE{k) =< *([5L*o > (2.63) 

Indeed, from the equation [dy — L)*o = and (2.53), it follows that 

dy < *o^*o >x^< %SL<i/o >^ . 

Taking the integral over y and using the following monodromy property of 5*o 

5^o(x, y + 1; /c) = W2ik){5^o(x, y; k) + 5E(k)^o). 
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we obtain (2.63). 

We begin by checking that the form < '^q5L A S'^q) >x is periodic in y. The shift of 
the argument y — > y + 1 gives 

< %6L A 5*0 >x — ^< %SL A 5*0 >a; + < %6L<ifQ >x A5E 

The second term on the right hand side can be rewritten as 5E A 5E and hence vanishes, 
due to the skew-symmetry of the wedge product. 

Next, we show that (b) is a consequence of the basic constraints defining the space 
C{b). Let *i be the formal Bfoch solution with the normalization *i(a;i,yi, k) = 1. Then 

*i(a;,y; k) = ^o{x,y;k)^^\xi,yi;k) 

and 

Resoo < *t5L A 5*1 > d/c = 2a; + Resoo ( SE A ^^oi^^^V^'^ 

V ^o{xi,yi;k) J 

In view of the constraints (2.30), we have 5E = 0{k~^). On the other hand, the second 
factor in the last term of the above right hand side also has order 0{k~^). The product 
has therefore order 0{k~'^) and its residue equals zero. 

To see that a; is a closed form, we express the operator L as 

n-2 
i=0 

which can be done in view of (2.57) and (2.58). Therefore 

2uj =< res9{D^-^6^^-^ A 5$) >= -S < resd{D^-'^6^) > . (2.65) 
and u) is closed. 

We turn now to the non-degeneracy of uj on C{b). Let V he a vector field such that 
u){Vi, V) — for all vector fields Vi. Let *i = 5*(F) be the evaluation of the one-form 
5* on V. Then the equality 

a;(Fi, V) = ReSoo < *o^i*i > d/c = 0, (2.66) 

holds for all degree n — 2 operators Li = dL{Vi). Since Li is arbitrary, it follows that 
*i = 0(A;~"^)*o- In view of (2.61) we have then 

6L{V)^o = {dy - L)*i = 0(A;-i)*o 
Hence 5L{V) = 0. This means that V = 0, and the non-degeneracy of oj is established. 
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It remains to exhibit the flows (2.41) as Hamiltonian flows. We recaU the classical 
definition of the Hamiltonian vector field dt corresponding to a Hamiltonian H and a two- 
form uj. The contraction i{dt)oj of uj with the vector field dt should be the one-form given 
by the differential of the Hamiltonian, i.e. the equality 

i{dt)u{X) = uj{X, dt) = dif(X), (2.67) 

should be fulfilled for all vector- fields X. 

The contraction of the form uj defined by (2.62) with the vector-field dm (2.41) is equal 

to 

2i{dm)(jo = KeSoo < %5Ldm'^ o > dk -ReSoo < %{dmL)5^o > dk. (2.68) 

Here we use the fact that the evaluations of the forms 5L and ^^^o on the vector field dm 
are equal by definition to 

SL{dm) = dmL, SMdm) = dm^O. 

From (2.44) it follows that 

Resoo < %SLdm'ifo > dk = ReSoo (< %SLAm^o > +^mik) < %SL'^o >) dk. 

The first term in the right hand side is zero due to the definition of ^'q- The usual formula 
for the implicit derivative 

6E{k)dk = -6k{K)dE, (2.69) 
implies that the second term is equal to 

-ReSoc^m{K)5k{K)dE = -ReSoo {K"^ + 0{K-^)) l^5Hsk-A dK"" = nSHn+m- 

\s=n J 

(2.70) 

(Recall that SHg — 0, s < n due to the constraints.) Consider now the second term in the 
right hand side of (2.68). The equation (2.41) for dmL and the defining equations for 
and imply 

ReSoo < *o(^m^)(^*o >x dk = Res^ody (< *^^^5*o >x)dk. 

Therefore, 

dy< *+yl^5*o >x dy = 5E{k) < *([Am*0 >x \y=Q- 

Jo 

The equality (2.37) implies 

<%Am^o >x= K^ + 0{K-^). 
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Hence, the second term in (2.68) is equal to 



Resoo < %{dmL)S'ifQ > dk = -nSHm+n- 

and Theorem 3 is proved. 

Example 1. For n = 2, the operator L is the second order differential operator of the 
form L = 9^ + u{x, y). The space Cipo) is the space of periodic functions with fixed mean 
value in x 

< u >x— Hi ^ < 5u >x— 0. 
and the symplectic form uj becomes 

1 r 

uj = —- < 5u A 5udx > 

This symplectic form reduces to the Gardner-Faddeev-Zakharov symplectic form when 
u{x, y) = u{x) is a function of a single variable x. In this case the KP equation reduces to 
the KdV equation. 

Example 2. For n = 3, the operator L is the third order differential operator L = 
d^ + udx +v. The space £(6o, 6i) is the space of doubly periodic functions u = u{x, y), v = 
v{x,y) satisfying the constraints 

< U >r— const, < V >r. — COUSt. 



The symplectic form u) works out to be 



1 / r 

UJ = — \ 6u A / 5vdx + Sv A 
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/ Sudx I 

'^Xo J 



In the case where u and v are functions of a single variable a;, this form gives a symplectic 
structure for the Boussinesq equation hierarchy 

_2 _2 

111 — ZVa; Uxx-i "^t — "^xx "^^xxx g'^'^a;- 

Note that the usual form of the Boussinesq equation, Utt + {^\uUx + \uxxx) ^ = 0, as an 
equation in one unknown function w, is the result of eliminating v from the above system. 

F. LgLX Equations 

In this section, we compare the results obtained in our formalism with the one- 
dimensional case, where the zero curvature equation reduces to the Lax equation, and 
where there is a rich theory of Hamiltonian structures. It turns out that the symplectic 
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structure constructed above reduces then to the so-called first (or generalized) Gardner- 
Faddeev-Zakharov symplectic structure. Thus our approach gives a new representation for 
this structure, as well as a new proof of its well-known properties (c.f. [10] [25]). As we shall 
see below, the second (Adler-Magri) symplectic structure requires a slight modification, 
which explains why it is special to the one dimensional case and has no analogue in the 
proposed Hamiltonian theory of two-dimensional systems. 

Our construction of the basic symplectic form uj easily extends to the construction of 
an infinite sequence of symplectic structures: 

Theorem 4. Let I he any integer > 0. Then the formula 

uj^i) = 1 Resoo E\k) < %6L A 5*0 > dk, (2.71) 

defines a dosed two-form on the space C{Hi, ■ ■ ■ , Hni-i) of doubly periodic operators L, 
subject to the constraints Hg = const, s = 1, . . . , n/ — 1. The equations (2.8) are Hamilto- 
nian with respect to this form, with the Hamiltonians (u) defined by (2.29). 

The proof of the theorem is identical to the proof for the basic structure u. Specializing 
to the subspace of periodic L with coefficients depending only on x, we can easily verify 
that the symplectic forms a;*-^-* coincide with the Gardner-Faddeev-Zakharov forms. 

The construction of the Adler-Magri Hamiltonian is less obvious, although formally 
it has the form (2.71) with I = —1 and the residue at infinity is replaced by the residue at 
E — 0. Let L be an ordinary linear differential operator of order n with periodic coefficients. 
Then for generic values of the complex number E, there exist n linearly independent Bloch 
solutions ^i{x, E) of the equation 

(L-E)*, = 0, (2.72) 

with different Bloch multipliers Wi{E), 

^i{x + 1, E) = Wi{E)^i{x, E). 

The value Pi{E) = log Wi{E) is called the quasi-momentum. Its differential dpi{E) is 
well-defined. (In our previous formal theory of Bloch solutions, there are also n different 
solutions corresponding to the same E, due to the relation E = ^ 0(A;"-2) = K"" which 
defines k and K only up to a root of unity.) We fix the Bloch solutions of the adjoint 
equation 

**(L-E) = 0. 

by the condition 

< ^*{x,E)^j{x,E) >a:=5ij 
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Theorem 5. The formula 

1 " 

'(-1) = - < ^*{x,0)6LASM^,0) >, (2.73) 



2 



where the constants Ri are given by 

Ri =< *|L(i)*i >, 

defines a closed two-form on the space of operators L with coefficients depending only on 
the variable x, and obeying the constraints Hg = const, s = 1, . . . ,n — 1. The equations 
(2.8) are Hamiltonian with respect to this symplectic form, with the Hamiltonians nHm{u) 
defined by (2.29). 

We observe that for n — 2, the operator L^^^ reduces to L^^-* = 2dx- The expression 
< W*L*^^)\E'i >x is then just the Wronskian of two solutions of the Schrodinger equation. 

The proof of Theorem 5 is analogous to the proof of Theorem 3. The formula (2.73) 
can be rewritten in the form 

u< >--- Reso — ^ ^^*^^^E)LW^,{x,E)>- ^^'^^^ 

Due to the summation over i, this expression is independent of the labeling of the Bloch 
functions. Thus on the right hand side, we have the residue of a well-defined function 
of E. The formula we need for the diff'erential of the branch of the quasi-momentum 
corresponding to the Bloch solution '^{x, E) of (2.72) is the following 

idp< <^*{x,E)L^^^-^i{x,E) >^dE. (2.75) 

Its proof is identical to the proof in the finite gap theory (see [36] ) . Consider the differential 
in the variable E of the Bloch function. Then 



(L - £;)d* = -*dE. 
Integrating from xq to xq + 1 the identity 

n 

= (**(L-E))d* = -(^**)d£; + J]a;^(**(LM^)). 

we obtain 

n 

dE = idp{^* [xo, E)L^'~>^,ixo, E) + Y, di-\^*ixo, i?)(L^*(xo, E))). 

J=2 
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The desired formula follows after averaging in xq this last identity. 

With the formula (2.75) for dE and the analyticity in E ior E ^ of all relevant 
expressions, we can, in the computation of the contracted form i{dm)u!^^\ reduce the 
residues at £■ = to the residue at E = oo and get the desired result. For example, we 
have 

-ResoJ2 < '^iSLdm-^i > = Resoo < ■^*{x, K)5Ldm-^{x, K) > ^ = nSHm- 

G. Higher Symplectic Structures in the Two-Dimensional Case 

In this section, we introduce higher Hamiltonian structures which exist in both one 
and two dimensions. We would like to emphasize that, in contrast with the previous 
results which have basically an algebraic nature, the following results require in general 
some additional assumptions on the long-time behavior of the solutions to the hierarchy 
of flows dmL = dtAm + [Am,L]. In these higher Hamiltonian structures, the role of the 
quasi- momentum k{K) for the basic structure is assumed by the quasi-momentums QmiK) 
corresponding to the higher times of the hierarchy. 

Our first step is to study in greater detail the quasi-momentums Q^niK), corresponding 
to higher times of the hierarchy. Their "densities" flrni^) made their first appearance in 
(2.44). They can be re-expressed as 

oo 

Qm{K) = Am'^oix, y, k)U=y=o = K"^ + J2 ^s,mK-' (2.76) 

The coefficients fls,m of Qm{k) are integro-differential polynomials in the coefficients of 
the operator L. As stated above, they do not depend on (x, y), but they do depend on the 
times t if the operator L evolves according to the equation (2.41), i.e. 0, = Q{K,t). Prom 
(2.44) and (2.45) it follows that 

diflm = dm^i. (2.77) 

Since the coefficients Ctg^rn ai'e independent of the choice of normalization point, they 
can be considered as functionals on the space of periodic operators L. 

The subsequent arguments are based on the variational formulas for these functionals, 
which were found originally in the case of finite-gap solutions in [36]. Following [36], we 
use the identity 

= {{%{L^^\5A^ + 5n^)*o) - {%{A^'^5L^q)) 

+ E (*o(^L'^^^''^ - L^')aH))5*o) , (2.78) 
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where L^-'^ and A^^^ are the descendants of the operators L and A defined by (2.60). Note 
that if L and A^ satisfy (2.44), then the equahty 

a^L(^) - dyA(£ + J^IL^'^), A^^-'^)] = (2.79) 

k 

holds. We now average (2.78) first in x and y, and then in the normahzation point Xq (the 
last averaging eliminates all terms with j > 1). The outcome is the equality 

SClmiK, t) < **L(^)*o >=< *o(^^'^'^^m - A^^^L)^o > +dm < *o^^'^5*o > • (2.80) 

Let , be the space of all periodic operators L which are stationary under the first 
mo-th fiow, i.e. which satisfy the condition 

dyAmo^[L,Am,]. (2.81) 

It should be emphasized that due to (2.77) onto this space the corresponding density of 
the corresponding quasi- momentum is a constant, i.e. does not depend on the times 

The space carrying a higher symplectic structure is a subspace V^l of , consisting 
of the stationary operators L satisfying in addition the following higher constraints 

oo 

n^^ = K-o + ^ n,^^^K-\ i}s,mo = Is, (2.82) 

s=l 

for a set / = (Ii, • • • , I^-i) of (n — 1) fixed constants. These constraints replace the 
constraints (2.30) of our previous considerations. The subspace V^l is invariant with 
respect to all the other flows corresponding to times U. 

Theorem 6. The formula 

= ^Resoo < m^i^l^L - L^^^SAmo) A S^o > dk, (2.83) 

defines a closed two- form on T>ml- The restrictions of the equations (2.41) to this space 
are Hamiltonian with respect to this form, with Hamiltonians nQrno,n+m- 

Remark 1. This statement has an obvious generalization if we replace the stationary 
condition (2.81) by the condition that L be stationary with respect to a linear combination 
of the first mo fiows, i.e. 

dyA=[L,Al A^^c,A,. 

i=0 
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Remark 2. For mo = 1, we have Ai = d, A^^ = 1, and the formula (2.83) becomes 
identical to (2.62), i.e. oji = oj. 

The proof of Theorem 6 is identical to that of Theorem 3, after replacing of (2.63) by 
the formula 

SQmoiK)dE =< %{L^^^SAmo - A^^lSL)^o > dk, (2.84) 

which is valid on Pmo- This formula is itself a direct corollary of (2.80) and (2.75). 

As an example, we consider the case n = 2, uiq = 3. For n = 2, the equations (2.41) define 
the KP hierarchy on the space of periodic functions u{x, y) of two variables. For A = L3, 
the condition (2.81) describes the stationary solutions of the original KP equation, i.e. the 
space of functions described by the equation 

Suyy + {6uux + Uxxx)x = 0. (2.85) 

Theorem 6 asserts that, besides of the basic Hamiltonian structure, the restriction of any 
flow of the KP-hierarchy to the space of functions u{x, y) subject to (2.85) is Hamiltonian 
with respect to the structure given by (2.83). 

In the one-dimensional case, the constraint (2.81) is equivalent to the restriction of the 
Lax hierarchy on the space of finite-gap solutions. This space is described by the following 
commutativity condition for the ordinary differential operators L and A of respective 
degrees n and mo 

[L,^^J = 0. (2.86) 

This condition is equivalent to a system of ordinary differential equations for the coefficients 
Ui{x) of L. Theorem 6 asserts that the restriction of the Lax hierarchy to the space of 
solutions to (2.86) is Hamiltonian with respect to the symplectic form (2.83). In particular, 
the first flow (which is just a shift in x) is Hamiltonian. For the KdV case the corresponding 
symplectic structure coincides with the stationary Hamiltonian structure found in [5] . 

Example 3. We return to the case n = 3 of Example 2, and consider this time operators 
A of order m = 2. Thus the operators L and A are given by L = d^ + ud^ + v^ A = d^ + ^u. 
The space T>c is the space of two quasi-periodic functions u{x,y) and v{x,y) satisfying the 
constraints 

< u >x= const, < V >x= const, < >x= const. 
The operators L^^^ and A^^^ are given by 

and the symplectic form ojrn of (2.83) becomes 

u =< -u5u A 5udx + 25v A / 5v dx + 25u A5v — -5ux A 5u > . 

^ J Xo J Xq ^ 
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H. Symplectic Structures under Ergodicity Assumptions 

It may be worthwhile to point out that the existence of the higher Hamiltonian struc- 
tures obtained in the previous section requires less that the stationary condition (2.81). 
The only item which was necessary to the argument was the possibility of dropping the 
last term (which was a full derivative in tmo) in the formula (2.80). 

This suggests considering the space 'D'f^l^ of all operators L with smooth periodic 
coefficients, for which the corresponding solution L{tjn^), L{0) = L of the equation (2.41) 
for m = niQ exists for all tmo with uniformly bounded coefficients. In this case we may 
introduce the quasi-momentum 

T 



1 r ~ 

0^o(K)= hm - / nrno{K,t)dtmo- 
T^oo 1 Jq 



Although in this definition, only the dependence on tmo eliminated through averaging, 
the quasi-momentum is actually also independent of all the other times ti, in view of 
(2.77). 

For L e the formula 

5nmo{K)dE =< %{L^^^SAmo - AglSL)^o >o dk, (2.87) 
where < f{x, y, t) >o stands for 

1 

< f{x,y,t)>o= hm - <f{x,y,t>dt, 

T^OO 1 Jq 

holds. Here we make use of the fact that if 6L is variation of the initial Cauchy data 1/(0) 
for (1.47) then the variation SL(tmo) is defined by the linearized equation 

dmo^L - dy6Amo + [SL, Amo] + [L, SAmo]- 

With (2.87), it is now easy to establish the following theorem 
Theorem 7. The formula 

ojmo = ^Resoo < niA^mlSL - L^^^SAmo) A <^*o >o dk, (2.88) 

defines a closed two-form on subspaces of T>^^ subject to the constraints (2.82). The 
restrictions of the equations (2.4.1) to this space are Hamiltonian with respect to this form, 
with the Hamiltonians nQ.mo,n+m- 

The space appears to be a complicated space, and we do not have at this moment 
an easier description for it. As noticed in [36], it contains (for an arbitrary mo) all the 
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finite-gap solutions. There exist a few other cases where we can justify the ergodicity 
assumption. For example, for the KdV hierarchy, the ergodicity assumption is fulfilled for 
smooth periodic functions with sufficiently rapidly decreasing Fourier coefficients. Indeed, 
if u{x) can be extended as an analytic function in a complex neighborhood of real values 
for X and y, 

\u{x)\<U, |Ima;|<g (2.89) 

then u{x, t) is bounded by the same constant for all t due to trace formulae. Using the ap- 
proximation theorem [37] for all periodic solutions to (2.5) (also called the KP-2 equation, 
by contrast with the KP-1 equation given by (3.19) below) by finite-gap solutions, we can 
prove the ergodicity assumption in the case when the Fourier coefficients Uij of u satisfy 
the condition \uij \ < C/g'*'"'"'-'!. 

Important Remark. In order to clarify the meaning of the higher symplectic forms and the 
higher Hamiltonians, it is instructive to explain its analogue for the usual Lax equations. 
The Lax equations dmL = [M^, L] obviously imply that the eigenvalues of L are integrals 
of motion, and usually they serve as Hamiltonians for the basic symplectic structure. The 
higher Hamiltonians correspond to the eigenvalues of the operator M^o instead of L. Of 
course, they are time dependent, but after averaging with respect to one of the times, 
namely t^^, they become nonlocal integrals of motion and can serve as Hamiltonians for 
the corresponding symplectic structures. 

I. The Matrix Case and the 2D Toda Lattice 

Our formalism extends without difficulty to a variety of more general settings. We 
shall discuss briefly the specific cases of matrix equations and of the Toda lattice, which 
correspond respectively to the cases where L is matrix-valued, and where the differential 
operator dx is replaced by a difference operator. 

Let L = ^^=QUi{x,y)dl. be then an operator with matrix coefficients Ui = (wf^) 
which are smooth and periodic functions of x and y, whose leading term tt"^ = u^SafS is 
diagonal with distinct diagonal elements ^ for a ^ P, and which satisfies w^^^i = 0. 
Then, arguing as in Section II. B, we can show that there exists a unique matrix formal 
solution = (*o^(a^) V'l k)) of the equation {dy — L)'^o = 0, which has the form 

oo n— 2 

y;k) = iI + J2 y)k-')eMkx + Unk'^y + Bi{y)k') (2.90) 

8=1 i=Q 

(where / is the identity matrix, = Hs'^) matrix functions, and Bi{y) — {B'^^{y)) = 
{Bf{y)dai3) are diagonal matrices), and which has the Bloch property 

<ifo{x + l,y;k) = •^o{x,y;k)wi{k), wi{k) = e^. 
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The formal solutfon ^o{x, y; k) has the Bfoch property with respect to y as well, 

*o(a;, y + l;k) = *o(a;, y; k)w2{k) 
with the Bloch multiplier W2{k) of the form 

CO n — 2 

W2{k) = (1 + 5^ J,/c-^)expK/c" + Yl ^^^')' 

s=l i=0 

where Jg and Bi are diagonal matrices. 

As noted at the end of Section II. B, the second Bloch multiplier defines the quasi- 
energy E{k). This defines in turn the functionals Si just as in (2.28), with the only 
difference the fact that they are now diagonal matrices. If we introduce the diagonal 
matrix K by the equality 

oo 

UnK'' = E{k) = log W2{k)=Unk''+ J] Sik'' (2.91) 

i=-n+2 

then we may define diagonal matrices Hg — (11^6^)3) by 

CO 

kI = K + J2 HsK-' (2.92) 

The definition of the commuting flows in the matrix case is then just the same as 
in the scalar case. The only difference is that the number of these flows is now times 
larger. The corresponding times are denoted hy t — (t^.m): and the flows are given by 

[dy - L, da,m - La,m] = (2.93) 

where L^t^^ is the unique operator of the form 

m — l 

La,m = Vadf + ^ Ui^(^a,m){t)dl, vf^ = 5a^(J^^ 
i=l 

which satisfies the condition 

i^a,m*o(a;, y; k) = ^o{x, y; k){vaK^ + 0{K-^)). (2.94) 

As before, the dual Bloch formal series ^'0(0;, /c) is defined as being of the form (2.58), 
and satisfying the equation 

ReSooTr(%vd^^o)dk = 0, m > (2.95) 
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for arbitrary matrices v. We have then 
Theorem 8. The formula 

u) = ^ReSooTr < %6L A > dk (2.96) 

defines a closed non- degenerate two-form on the space of periodic operators L subject to 
the constraints = constant, (3 = 1, ■ ■ ■ , s = 1, ■ ■ • , n — 1. The equations (2.93) are 
Hamiltonian with respect to this form, with Hamiltonians nH^_^_^ defined by (2.92). 

Example 4- Consider the case n = 1, where the operator L is of the form L = 
Adx + u{x, y), with A the N x N matrix A"'^ = aoi?>oL,f3 and u{x, y) is an N x N matrix 
with zero diagonal entries = 0. In this case the symplectic form (2.96) becomes 



(jj 



= V A 6u'''f^. 



aa — ap 



Finally, as a basic example of a system corresponding to an auxiliary linear equation 
where the differential operator dx is replaced by a difference operator acting on spaces of 
infinite sequences, we consider the 2D Toda lattice. 

The 2D Toda lattice is the system of equations for the unknown functions (fn = 

It is equivalent to the compatibility conditions for the following auxiliary linear problem 

We consider solutions of this system which arc periodic in the variables n and y = 
(t_i_ + t-). The relevant linear operator is the difference operator 

with periodic coefficients Vn{y) = Vn+N{y) = Vn{y + I) and Cn{y) = Cn+N{y) = Cn{y + 1). 
Then, arguing as in Section II. B we can show that there exist unique formal solutions 
= *i'^''(2/; k) of the equation 

{dy - L)*(±) = 0, 

which have the form 

¥,^\y;k) = (^6^\n,y)k-''^ ^^v+Biv) ^ ^(+) ^ ^2.97) 
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the Bloch property 

and which are normahzed by the condition 

The coefficients ^i^-* can be found recursively. The initial value ^q^^ = 1 and the condition 
that ^['^^ is periodic in n define the function B{y) in (2.97) 



The only difference with the previous differential case is the definition of the leading term 
\n,y). Let us introduce (Pn{y) — log^g \''^iy)- Then we have 

The periodicity condition for -* requires the equality 

N-l 

^ dy^Pn = 0, 
n=0 

which allows us to define (fin uniquely through Cn' 

n-l / N-l 



(Pn = ^^ogCn- I ( ^ ^ logCn ) dy 

n=l V n=0 / 



The formal solutions ^^\y,k) have the Bloch property with respect to y as well, 

^i^\y+l;k) = ^i^\y;k)wi'^\k) 
with the Bloch multipliers W2'^\k) of the form 



vs=0 



As noted at the end of Section II. B, the second Bloch multiplier defines the quasi-energy 
E^^\k) and the functionals ^j^^"* just as in (2.28). If we introduce the variable 

CO 

K = E^^\k) = \ogw^^\k) ^k + ^el'^h-' 

i=0 
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then we may define the functionals i^i^^ by 

logA; = logK± (2.98) 

s=0 

The definition of the commuting flows in the discrete case is then the same as in the scalar 
case. The basic constraints that specify the space of periodic functions Vn and Cn have the 
form 

n N 

Jo = ^ log Cn(y) = const, h = ^2 '"niy) = const. (2.99) 

n=l n=l 

The corresponding times are denoted hy t = (t±,m) and the flows are given by 

[dy-L,d±,m-L±,m]=0, (2.100) 
where L±^rn is the unique operator of the form 

m 
i=0 

which satisfies the condition 

L±,^vE'i±) = vE'i±)(i^"^ + 0(i^-i)). 
The dual formal series are defined as the formal series of the form 

which satisfy the equations 

for all integers m. 
Theorem 9. The formula 

dk 

u = Resoo < A 5*^+) - **'(-)5L A S^^'^ > — =< {25vn - Sdy(pn) A S(fn > 

(2.101) 

defines a closed non-degenerate two- form on the space of periodic operators L subject to 
the constraints (2.99). The equations (2.100) are Hamiltonian with respect to this form, 
with Hamiltonians defined by (2.98). 
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Finally, we would like to conclude this Chapter by calling the reader's attention to 
[59] [60], where symplectic structures are discussed from a group theoretic viewpoint, and 
where applications of soliton theory to harmonic maps are given. 



III. GEOMETRIC THEORY OF 2D SOLITONS 



In Chapter I, we had developed a general Hamiltonian theory of 2D solitons. The 
central notion was the symplectic form (1.1), which was defined on the infinite-dimensional 
space C{h) of doubly periodic operators obeying suitable constraints. Our main goal in 
this Chapter is to present and extend the results of [39]. In this work, as described in 
the Introduction, a natural symplectic form lom was constructed on Jacobian fibrations 
over the leaves of moduli spaces A4 g{n,m) of finite-gap solutions to soliton equations. 
Imbedded in the space of doubly periodic functions, the form was shown to coincide 
with uj (this was in fact our motivation for constructing a general Hamiltonian theory 
based on iv in this paper). Although the infinite-dimensional symplectic form iv and its 
variants in Chapter II can be expected to play an important role in an analytic theory of 
solitons, it is the geometric and finite-dimensional form ivm which has provided a unifying 
theme with topological and supersymmetric field theories. 

In Chapter II, we have seen how a differential operator L determined a Bloch function 
^^0) which was a formal series in a spectral parameter k or K. The key to the construction 
of finite-gap solutions of soliton equations is the reverse process, namely the association of 
an operator L to a series of the form of ^'o- To allow for evolutions in an arbirary time 
trm it is convenient for us to incorporate a factor e^'"'^ in for each t^, and consider 
series ^'o(^; k) of the form 



As usual, all the times U except for a finite number have been set to 0. Then the operators 
Lm are uniquely defined by the requirement that 



(this is equivalent to the earlier requirement that (L^ — k'^)^Q = O{k~^)^o in the case 
of ^s{t) independent of tm)- In particular, we have the following identity between formal 
power series 




(3.1) 



(5, 



'm 



CO 




(3.2) 
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This identity assumes its full value when the formal series ^o{t; k) is a genuine convergent 
function of k and has an analytic continuation as a meromorphic function with g poles on 
a Riemann surface of genus g. In this case, the equation (3.2) with zero right hand side 
becomes exact. The null space of [dn — dm — Lm\ is parametrized then by k and is 
infinite-dimensional. Since [dn — Ln, dm — Lm] is an ordinary differential operator, it must 
vanish. Thus a convergent \I/o(t; k) gives rise to a solution of the zero curvature equation 
[dn — Ln, dm — Lm] = 0. The algebraic-geometric theory of solitons provides precisely the 
geometric data which leads to convergent Bloch functions. These functions are now known 
as Baker-Akhiczer functions. 

A. Geometric Data and Baker- Akhiezer Functions 

In a Baker-Akhiezer function, the spectral parameter k is interpreted as the inverse 
k = of a local coordinate z on a Riemann surface. Thus let a "geometric data" 
(r, P, z) consist of a Riemann surface F of some fixed genus (7, a puncture P on F, and a 
local coordinate k~^ near P. Let 71, ■ ■ ■ , 7^ be (7 points of F in general position. Then for 
any t = (ti)^i, only a finite number of which are non-zero, there exists a unique function 

z) satisfying 

(i) \I/ is a meromorphic on F \ P, with at most simple poles at 7i, ■ ■ ■ , 7^; 

(ii) in a neighborhood of P, \1/ can be expressed as a convergent series in k of the form 
appearing on the right hand side of (3.1). 

The exponential factor in (3.1) describes the essential singularity of ^{t;z) near P. 
Alternatively, we can view it as a transition function (on the overlap betwwen F \ P and 
a neighborhood of P) for a line bundle C{t) on F. The Baker-Akhiezer function is then 
a section of >C(t), meromorphic on the whole of F. 

The form of the essential singularity of implies that has as many zeroes as it 
has poles (equivalently, the line bundle jC{t) has vanishing Chern class). Indeed, d'^/'^ = 
d{J2°^i tik~'^) + regular, and thus has no residue at P. From this, the uniqueness of 
the Baker-Akhiezer function follows, since the ratio ^ of two Baker-Akhiezer functions 
would be a meromorphic function on the whole of F, with at most g poles (corresponding 
to the zeroes of \E'). By the Ricmann-Roch Theorem, it must be constant. Finally, the 
existence of \1/ can be deduced most readily from an explicit formula. Let A^, - ■ ■ , Ag 
Bi, ■ ■ ■ , Bg he a canonical homology basis for F 



(a) dujj, Tjk be respectively the dual basis of holomorphic abelian differentials and the 
period matrix 



Aj nAk = Bj nBk = 0, Aj nBk = Sjk, 



(3.3) 



and let 
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(b) 0{z\t) the Riemann ^-function; 

(c) df2° the Abehan differential of the second kind with unique pole of the form 

dn^i = {l + 0{k-'-^))dk' (3.4) 
normalized to have vanishing A-periods 

df]0 = 0, 1 < / < ^; (3.5) 



/ 

J A 



(d) Po a fixed reference point, with which we can define the Abel map A: z eT ^ A{z) e 

and Abelian integrals by 

Ajiz)= Tdo;,-, n'^iz)= TdJ^o (3.6) 

JPo JPo 

Here the Abel map as well as the Abelian integrals are path dependent, and we need 
to keep track of the path, which is taken to be the same in both cases. 

(e) Z = K — X]f=i ^(Ts)) where K is the vector of Riemann constants (c.f. [48]). 

Using the transformation laws for the ^-function, it is then easily verified that the 
following expression is well defined, and must coincide with '^{t; z) 

In this formalism, the role of the quasi-momentum p is now assumed by the Abelian 
integral and we just define p = With the Baker- Akhiezer function assuming now 
the former role of formal Bloch functions, we can construct a hierarchy of operators 
as in (2.37). The requirement that {dm — Lm)^{t'i z) = 0{z)'^{t\ z) determines recursively 
the coefficients of Lm as differential polynomials in the ^g. The crucial improvement over 
formal Bloch functions is that here, this requirement actually implies that 

(a„-L„)*(t;^) = (3.8) 

identically. In fact, [dn — Ln)'^{t;z) satisfies all the conditions for a Baker- Akhiezer 
function, except for the fact that the Taylor expansion of its coefficient in front of the 
essential singularity exp(^^-^ t^/c') starts with k~^. If this function is not identically 
zero, it can be used to generate distinct Baker-Akhiezer functions from any given one, 
contradicting the uniqueness of Baker-Akhiezer functions. As noted earlier, (3.8) im- 
plies that [dn — Ln,dm — Lm]^{t',z) = 0, and hence that the zero curvature equation 
[dn — Ln, dm — Lm] = holds. In summary, we have defined in this way a "geometric map' 
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Q which sends the geometric data (F, P, 2;; 71, • • • , 7^) to an infinite hierarchy of operators 
[33] [34] 

6;:(r,P,^;7i,..-,7,)^(Lj-2 (3.9) 

The expression (3.7) for \l/(t;z) leads immediately to explicit solutions for a whole 
hierarchy of soliton equations. Let ti = x, t2 = y, = t, and n = 2. We find then 
solutions u(x, t) of the KP equation (2.4) expressed as [34] 

u{x,y,t) ^2dllog 9 {x f dQ^ + y f dQ^ + t f dft^ + Z\t) + const. (3.10) 

Jb J b J b 

This formula is at the origin of a remarkable application of the theory of non-linear inte- 
grable models, namely to a solution of the famous Riemann-Schottky problem. Accord- 
ing to the Torelli theorem, the period matrix defines uniquely the algebraic curve. The 
Riemann-Schottky problem is to describe the symmetric matrices with positive imaginary 
part which are period matrices of algebraic curves. Novikov conjectured that the function 
u{x, y, t) = 2dl\og9{Ux + Vy + VFt|r) is a solution of the KP equation if and only if the 
matrix r is the period matrix of an algebraic curve, and U, V, W are the i?-periods of the 
corresponding normalized meromorphic differentials with poles only at a fixed point of the 
curve. This conjecture was proved in [3] [55] . 

The dual Baker- Akhiezer function 

For later use, we also recall here the main properties of the dual Baker- Akhiezer func- 
tion \I/+(t; z) which coincides with the formal dual series defined in Section II. D. To define 
"^^{t; z), we note that, given g points 71, ■ ■ ■ , 7g in general position, the unique meromor- 
phic differential dO = d{z~^ + Y1^2 '^sZ^) with double pole at P and zeroes at 71, ■ ■ ■ , 7^, 
must also have g other zeroes, by the Riemann-Roch theorem. Let these additional zeroes 
be denoted by 7^, ■ ■ ■ , 7^. Then the dual Baker- Akhiezer function '^^{t] z) is the unique 
function \I/^(t; z) which is meromorphic everywhere except at P, has at most simple poles 
at 7]^, • • • , 7+, and admits the following expansion near P 

^+{t;z) = exp{-J2tnZ-n{l + J2Cit)zn 

n=l s=l 

To compare the dual Baker- Akhiezer function with the formal dual Bloch function ^'q' of 
Section II. D, it suffices to observe that 

Resp^+{t-z){d^-^{t-z))dn = 0, 

since the differential on the left hand side is meromorphic everywhere, and holomorphic 
away from P. Together with the normalization ^^+(0; z) = 1, this implies that '^'^ indeed 



36 



coincides with the formal dual function '^q. An exact formula for '^'^{t; z) can be obtained 
from (3.7) by changing signs for t and by replacing the vector Z by Z'^. From the definition 
of the dual set of zeroes 7^*", • • • , 7^", this vector satisfies the equation Z + Z^ = 2P + K, 
where K is the canonical class. Recalling that the quasi- momentum p was defined to be 
p = ill, we also obtain the following formula for the differential dVL we introduced earlier 

dn ^ ^ . (3.11) 

The Multi-Puncture Case 

The above formalism extends easily to the case of N punctures Pq, (with one marked 
puncture Pi). The Baker- Akhiezer function is required then to have the essential sin- 
gularity 

oo oo 
s=0 i=l 

where are local coordinates near each puncture Pq,, tia are given "times", only a finite 
number of which are non-zero, and the coefficient at Pi is normalized to be 1 for s = 0. 
We can introduce as before dO,^^ associated now to each puncture Pq, and their Abelian 
integrals fl^^. Then the Baker- Akhiezer function ^{t;z) becomes 

. JiMz) + 2^ Ell §B + z\T)e{A{p,) + z) 
' e{A{z) + z)e{A{p,) + ^ Eq=i TZi §b '^^L + z\t) 

^ N oo N 

xexp J2J2^^-^ic.izn (3-13) 

^a=li=l ^ 

For each pair (a, n) there is now a unique operator Lan of the form (2.36) so that (dan — 
Lan)i^it, z) = 0, with dan = d/dtan- The operators Lan satisfy the compatibility condition 

[dan ~ -t'an) ^/3m ~ -^'/Jm] = 0- 

Periodic Solutions 

In general, the finite-gap solutions of soliton equations obtained by the above construc- 
tion are meromorphic, quasi-periodic functions in each of the variables tai (a quasi-periodic 
function of one variable is the restriction to a line of a periodic function of several vari- 
ables). We would like to single out the geometric data which leads to periodic solutions. 
For this we need the following slightly different formula for the Baker- Akhiezer function. 

Let dQia be the unique differential with pole of the form (3.4) near Pq,, but normalized 
so that all its periods be purely imaginary, and define the function $((^1, • • • , (2g', z) by the 
formula 

*(C; z) = /j/ 7\nrAfD\\r TT ^ . exp(27r^ > Ak{z)Ck+g) 3.14 

0{A{z) + Z)0{A{Pi) + Ckek + Ck+gTk + Z\t) ^ 
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where = (0, • • • , 0, 1, 0, • • • , 0) are the basis vectors in C^', and Tk are the vectors with 
components Tjk- We observe that $ is periodic of period 1 in each of the variables 
Ci) ■ • ■ ) C23- Then the Baker- Akhiezer function can be expressed as 



(3.15) 



where we have denoted by Uia the real, 25f-vector of periods of dfij, 




2TTi 



1 




In particular, for geometric data {F, Pa-, z^} satisfying the condition 



1 




(3.16) 



the Baker- Akhiezer function is a Bloch function with respect to the variable t if we set tia = 
aicft, with Bloch multiplier w = exp(^^^ aia^iaiz)). The coefficients of the operators Lo-n 
are then periodic functions of t. As an example, we consider the one-puncture case. If we 
express the data under the form Ui — 27rmfc//i, C/^ — 27rnfc/Z2, with mk,nk G Z, then the 
corresponding solution of the KP hierarchy is periodic in the variables x = ti, y — t2, with 
periods h and I2 respectively. 

Real and Smooth Solutions 

There are two types of conditions which guarantee that the solutions obtained by the 
above geometric construction are real and smooth for real values of tia- We present them 
in the case of the KP hierarchy. 

Assume that the geometric data defining the Baker- Akhiezer function is real, in the 
sense that 

(a) the algebraic curve F admits an anti-holomorphic involution i : F — > F; 

(b) the puncture Pi is a fixed point of t; 

(c) the local coordinate in a neighborhood of Pi satisfies the condition k{L{z)) = k{z); 

(d) the divisor (71, • • • ,7^) is invariant under l, i.e., i{'js) = 7ct(s)) where u is a permuta- 



This is an immediate consequence of the uniqueness of the Baker-Akhiezer function and 
the fact that both sides of the equation have the same analytic properties. In particular, 
the coefficients of L„ and the corresponding solutions of the KP hierarchy are real. 



tion. 

Then the Baker-Akhiezer function satisfies the reality condition 



^{t;i{z))^^{t;z) 



(3.17) 
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In order to have real and smooth solutions, it is necessary to restrict further the 
geometric data. In general, the set of fixed points of any anti-holomorphic involution on 
a smooth Riemann surface is a union of disjoint cycles. The number of these cycles is less 
or equal to g + 1. The algebraic curves which admit an anti-involution with exactly g + 1 
fixed cycles arc called M- curves. We claim that the coefficients of Ln are real and smooth 
functions of all variables ti when F is an M-curve with fixed cycles Aq, Ai, ■ ■ ■ , Ag, and 
P G Aq, 7s G As, s = 1, ■ ■ ■ , g. To see this, we note that, from the explicit expression for 
the Baker- Akhiezer function, the coefficients of L„ have poles at some value of ti if and 
only if 

e{A{Pi) + J2 Uiti + Z) = (3.18) 

i 

The monodromy properties of the ^-function imply that the zeros of the function 6{A{z) + 
Si Uiti + Z) are well-defined on F, even though the function itself is multi-valued. The 
number of these zeroes is g. They coincide with the zeroes of 'i!{t,z). In view of (3.17), 
the Baker- Akhiezer function is real on the cycles Ag. On each of the cycles Ai, • • ■ , Ag, 
there is one pole of There must then be at least one zero on the same cycle. Hence 
all zeroes of 0{A{z) + Uiti + Z) are located on cycles Ag. Since Pi G the equation 
(3.18) cannot be fulfilled for real values of ti. 

We observe that the real and smooth solutions of the KP hierarchy corresponding to 
M-curves with a fixed puncture, are parametrized by the points of a real ^-dimensional 
torus which is the product of the g cycles Ag. If we choose these cycles (as our notation 
suggests) as half of a canonical basis of cycles, then this torus is the real part of the 
Jacobian. 

In the theory of real and smooth solutions the equation (2.4) is called the KP-2 
equation. The other, so-called, KP-1 equation is the other real form of the same equation. 
It can be obtained from (2.4) by changing y to iy. Thus the KP-1 equation is given 
explicitly by 

~'^^yy ~ i^t ~ 2^^^ ~ '^'^xxx)x- (3.19) 

As complex equations (2.4) and (3.19) are equivalent. But the conditions which single out 
real and smooth solutions are different. These conditions for the KP-1 equation may be 
found in [37] . Brieffy they are : 

Assume that the geometric data deffning the Baker- Akhiezer function is real, in the 
sense that 

(a) the algebraic curve F admits an anti-holomorphic involution i : F — > F; 

(b) the puncture Pi is a ffxed point of t; 

(c) the local coordinate in a neighborhood of Pi satisffes the condition k{L{z)) = 
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(d) the divisor (71, • • • , 7^) under l becomes the dual divisor 7i", • " " ) 7^ i-^-) i^ila) = 7j"( 
where cr is a permutation. 

Then the Baker- Akhiezer function satisfies the reahty condition 



(3.20) 



where the new variables t' = {t'l, . . .) are equal to t2m+i — ^2m+i, ^2m = ^^2m- As before, 
this is an immediate consequence of the uniqueness of the Baker- Akhiezer function and 
the fact that both sides of the equation have the same analytic properties. In particular, 
the coefiicients of L„ and the corresponding solutions of the KP-1 hierarchy are real for 
real values of t'. 

The further restriction of geometric data corresponding to real and smooth solutions 
of the KP-1 hierarchy is as follows. The fixed cycles ai, • • • ,a/ of t should divide F into 
two disconnected domaines r='=. The complex domain defines the orientation on the 
cycles considered as its boundary. The differential dfl of (3.11) should be positive on 
with respect to this orientation. 

B. Moduli Spaces of Surfaces and Abelian Integrals 

The space {F, P, z, 71, ■ ■ ■ , 7^} provides geometric data for solutions of a complete 
hierarchy of soliton equations, and is infinite-dimensional. In the remaining part of this 
paper, we concentrate rather on a single equation of zero curvature form [dy — L,dt — 
A] = 0. The geometric data associated with the pair (L, A) corresponds to the Jacobian 
bundle over a finite-dimensional moduli space A4g{n,m) of Riemann surfaces with a pair 
of Abelian integrals {E, Q) with poles of order n and m respectively at the puncture P. 
The associated operators (L, A) are then operators of order n and m, and are obtained by 
the basic construction (3.9), after imbedding M.g{n, m) in the space (F, P, z) of geometric 
data. Alternatively, we may choose to represent the equation [dy — L,dt — A] = as a 
dynamical system on a space of operators L, with t as time variable. In this case, a finite- 
dimensional and geometric space of operators L is obtained by the same construction as 
just outlined, starting instead from the Jacobian bundle over the moduli space Ai g{n) of 
Riemann surfaces F with just one Abelian integral with pole of order n at the puncture. 

More precisely, given (F, E), a geometric data (F, P, z) is obtained by setting the local 
coordinate z = near the puncture P to be 



where n > 1 and R are respectively the order of the pole of E and its residue at P. When 
n = 0, we set instead 



£; = ^-^ + i?^log^ 



(3.21) 



E = R'^logz. 



(3.22) 
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This gives immediately a map 



{r,E)^{T,P,z) 
(r,£;,7i,---,7j^(r,P,2;,7i,---,73)^ L (3.23) 

where the operator L is characterized by the condition {dy — L)\l/ = 0, with \E'(a;,y;/c) 
the Baker- Akhiezer function having the essential singularity exp(A;x + k^y), k = z~^. In 
presence of a second Abelian integral Q, we can select a second time t, by writing the 
singular part Q+{k) of Q as a polynomial in k and setting 

Q+{k) = aik+ ha^/c"" 

ti = ait, l<i<m. (3.24) 

This means that we consider the Baker-Akhiezer function '^{x,y,t;k) with the essential 
singularity exp{kx + k'^y + Q +{k)t) , and construct the operators L and A by requiring that 
{dy — L)\E' = {dt — A)^ = 0. The pair (L, A) provides then a solution of the zero-curvature 
equation. By rescaling t, we can assume that A is monic. Altogether, we have restricted 
the geometric map Q of (3.9) to a map on finite-dimensional spaces, which we still denote 

6;:(r,£;;7i,---7,)->(r,p,^)^ (l) 

g:{r,E,Q;-fi,---fg)^{r,P,z,t)^ {L,A) (3.25) 

Here we have indicated explicitly the choice of time in the geometric data. The proper 
interpretation of the full geometric data (F, Q; 71, • • -7^) is as a point in the bundle 
M^{n,m) over M.g{n,m), whose fiber is the ^-th symmetric power S^{T) of the curve. 
The ^-th symmetric power can be identified with the Jacobian of T via the Abel map 

(71, • • • , 75) ^ 0j = XI / ^^j' (^-2^) 

More generally, we can construct the bundles Ng{n,m) and Mgin) with fiber S^{V) over 
the bases A4g{n, m) and A4g{n) respectively 

Sk{T) U^{n,m) S^{T) A/J(n) 

i i (3.27) 

Mg{n,m) Mg{n) 

Thus the bundles J\fg^^{n, m) = Mg{n, m), and Mg^^{n) = Mg{n) are the analogues in the 
our context of the universal curve. Returning to soliton equations, the geometric map Q of 
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(3.25) can now be succinctly described as a map from the fibrations M^{n) and m) 
into the spaces respectively of operators L and pairs (L, A) of operators 



g:J\f^{n,m)^{L,A) (3.28) 

We emphasize that, although the operators in its image are not all periodic operators, the 
ones arising later upon restriction of Q to suitable subvarieties of J\f^{n,m) and M^{n) 
with integral periods (c.f. Section III.C) will be. 

We conclude this section by observing that, in the preceding construction, L and 
A depend only the singular part of Q, and hence are unaffected if dQ is shifted by a 
holomorphic differential. As we shall see below, the appropriate normalization in soliton 
theory is the real normalization by which we require that 

Re ^ dQ = (3.29) 

In the study of N=2 supersymmetric gauge theories, holomorphicity is a prime considera- 
tion, and we shall rather adopt in this context the complex normalization 

f dQ ^ 0, 1 < j < g (3.30) 

We note that each normalization provides an imbedding of A4g{n) into A4g{n, 1), by 
making the choice Q+{K) — K, with periods satisfying either (3.29) or (3.30), so that the 
operator A is just A = either case. The image of M.g{n) in M.g{n, 1) does depend 

on the normalization, however. 

C. Geometric Symplectic Structures 

We begin by discussing the basic local geometry of the moduli spaces Aig{n) and 
A4g{n,m). They are complex manifolds with only orbifold singularities, of dimensions 

N 

diuiMgin) = - 3 + 2N + y^^rig 

a=l 

TV 

dim Mg{n,m) = 5g - 3 + 3N + ^{ua + rua) (3.31) 

Indeed, the number of degrees of freedom of an Abelian integral E with poles of order 
n = {ua) is 1 + J2a=i(''^a + 1) — l + g = N + g + J2a=i "^oii where the first 1 corresponds 
to the additive constant, and the remaining integer on the left is the dimension of mero- 
morphic differentials with poles of order < tt-q, + 1 at each P^. For ^ > 1, the dimension 
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of the moduli space of Riemann surfaces with N punctures is 3g — 3 + which leads 
immediately to (3.31). For ^ < 1, it is easily verified that the same formula (3.31) holds, 
although the counting has to incorporate holomorphic vector fields and is slightly different 
in intermediate stages. 

We can introduce explicit local coordinates on M.g{n,m). To obtain well-defined 
branches of Abelian integrals, we cut apart the Riemann surface T along a canonical 
homology basis Ai,Bj, i,j = l,---,g, and along cuts from Pi to Pa for each 2 < a < 
N. Locally on M.g{n,m), this construction can be carried out continuously, with paths 
homotopic by deformations not crossing any of the poles. Denote the resulting surface by 
Tcut- On Tcitt, the Abelian integrals E and Q become single- valued holomorphic functions, 
and we can introduce the one-form dA by 

dA = QdE (3.32) 

We observe that dA has a singularity of order Ua + rua -|- 1 at each puncture Pa- Now the 
Abelian integral E defines a coordinate system Za near each Pa by 

E = z-""- +R^\ogZa (3.33) 

when Ha is strictly positive. When na = 0, we write instead 

E = R^\ogZa. (3.34) 

The coordinate Za can be used to fix the additive normalization of the Abelian integral 
A, and to describe its Laurent expansion near each puncture. Thus we fix the additive 
constant in A by demanding that its expansion in zi near Pi have no constant term. The 
parameters Ta,ii 1 < i < ria -\- ma-, -Rq, 2 < a < N , can then be defined by 

Ta,i = --Resp^ (^adA), 1 < a < iV, 1 <i <na + ma 

i?^ = RespJdA) (3.35) 

The parameters Ta,i (1 < a < N) and Ra (2 < a < N) account for Y^^^iina+ma) + N-l 
parameters. The remaining parameters needed to parametrize A4g{n,m) consist of the 
- 2 residues of dE and dQ 

i?f = Resp„d£;, r2 = Resp„dQ, a = 2,---,N (3.36) 

and the following 5g parameters which account for the presence of non-trivial topology 



TAi,E = f dE, TB,,E = <h dE (3.37) 



TA„Q = f dQ, TB„Q = f dQ (3.38) 
JAi JBi 

(p QdE, i = l,---,g (3.39) 



Ai 
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Theorem 10. Let T> he the open set in M.g{m,n) where the zero divisors of dE and dQ, 
namely the sets {z;dE{z) = 0} and {z;dQ{z) = 0}, do not intersect. Then 

(a) Near each point in T>, the 5^ — 3 + 3N + X]a=i('^a + "^a) parameters R^, R^, R^, 

Ta,k, '^Ai,E, T~Bi,E, TAi,Q, TBi,Q, di have linearly independent differentials, and thus 
define a local holomorphic coordinate system for Aigin^m) ; 

(b) The joint level sets of the set of all parameters except define a smooth g-dimensional 
foliation of T>, independent of the choices we made to define the coordinates themselves. 

Theorem 10 is proved in [39]. Since M.g{n) can be imbedded in Aig{n, m) by choosing 
Q to have Laurent expansion Q+{K) ~ K and fixing its periods, Theorem 1 also 
provides local coordinates for Mg{n). Specifically, the coordinates of Mg{n) which arise 
this way are 

• Tq, j with 1 < i < ria for a > 2, Ti j with 1 < z < ?ii — 1 (since the normalization 
Q — + 0{K) fixes the coefficients of the two leading terms Ti and Ti ^^ in 
the singularity expansion of A near Pi to be Ti = r^Xi ' ~ 

• the 2N — 2 residues of dE and dA at for a > 2; 

• the Bk periods of dQ, the Ak and periods of d-E, and the A^ periods of dA, 
for a total of - 3 + 2iV + Y.a=i ^a- 

The intrinsic foliation obtained in Theorem 10 is central to our considerations, and we 
shall refer to it as the canonical foliation. Our goal is to construct now a symplectic form 
u! on the complex 2^-dimensional space obtained by restricting the fibration Af^{n, m) to a 
^-dimensional leaf M. of the canonical foliation of Alg(n, m). For this, we need to extend 
the differentials dE and dQ to one-forms on the whole ffbration, and distinguish between 
the two ways this can be done. 

One way is to consider the Abclian integrals E and Q as multi-valued functions on the 
ffbration. Despite their multivaluedness, their differentials along any leaf of the canonical 
ffbration are well-deffned. In fact, E and Q are well-deffned in a small neighborhood of the 
puncture Pi. The ambiguities in their values anywhere on each Riemann surface consist 
only of integer combinations of their residues or periods along closed cycles. Thus they 
are constant along any leaf of the canonical foliation, and disappear upon differentiation. 
The differentials along the ffbrations obtained this way will be denoted hy 5E and 5Q. 
Restricted to vectors tangent to the fiber, they reduce of course to the differentials dE and 
dQ. 

The other way is to trivialize the fibration by using an Abelian integral, say as local 
coordinate in F (c.f. (3.21) and (3.22)). Equivalently, we note that at any point of jV", the 
varieties E = constant are intrinsic and transversal to the fiber. Thus any one-form dQ 
on the fiber can be extended to a one-form on the total space of the fibration, by making 
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it zero on vectors tangent to these varieties. We still denote this one-form by dQ. In the 
case of the Abelian integral E, we have 5E = dE, but this is not true in general. To 
compare dQ with SQ, let ai, • • • , be local coordinates for the leaf. Then ai, ■ ■ ■ , ag, E 
are local coordinates for the whole fibration, dai,- • - jdag, dE are a basis of one-forms, 
and dQ = j^dE. On the other hand, 

5Q = dQ + y ^dai = dQ + 5^Q. (3.40) 

Similarly, the full differential 6{QdE) on the total space of the fibration is well-defined, 
despite the multivaluedness of Q. The partial derivatives ^^^(Qdi?) along the base are 
all holomorphic, since both the singularities and the ambiguities in the differential are 
constant, and disappear upon differentiation. Recalling that da;^ denotes the basis of 
holomorphic differentials dual to the homology basis Ak, Bk, we can then write 

■^{QdE) = du;i. (3.41) 

By extending this construction to the fibration of Jacobians over a leaf M. of the canonical 
foliation, we obtain the desired geometric symplectic form. Furthermore, this symplectic 
form coincides with the symplectic form constructed in Section II. E, upon imbedding the 
fibration of Jacobians in the space of soliton solutions of the equation [dy — L,dt — A] = 
by the geometric map Q of (3.28) [39]: 

Theorem 11. (a) The following two-form on the fibration J\f^{n,m) restricted to a leaf 
Ai of the canonical foliation of A4g{n,m) 

9 9 9 

OJM = 5(5]Q(7i)d£;(7i)) = Y,^Q{li) ^dE{-i,) = J]dai Ada;^ (3.42) 

i=l i=l i=l 

defines a symplectic form; (b) Under the geometric correspondence Q, we have 

U}M = Q*{<^m) (3.43) 

where u>m is the symplectic form constructed in Theorem 7. 

D. Generalizations and Extensions 

The above set-up can be easily adapted to a variety of important equations. 

For the 2D Toda lattice, the differential operator dx is replaced by a difference op- 
erator. Thus for the geometric data, we consider the moduli space of Riemann surfaces 
F with two punctures P±, and Abelian differentials E, Q which are real-normalized, with 



45 



E having a pole of order 1 at P+, while dQ has simple poles at P±. This is a leaf in the 
foliation of M.g{l, 0; 0, 0). The discrete analogue of the Baker-Akhiezer function '^{z; t) is 
now given by a sequence ^riit+, t-] z) of Baker-Akhiezer functions characterized by 

oo 

^nit+,t-;z) = z^''{J2Un,t)zi)eMt+z+^ +t-zZ^) (3.44) 

As in Section II. I, we define the operators L and A as difference operators acting on 
sequences ipn which satisfy {dt_^_ — L)'^ = {dt_ — A)'^ — 0. Their coefficients c^, Vn, (j)n are 
now functions of the geometric data. In analogy with (3.43), we have then 

^G*{i^M) = - X] ^^^P± < A > dp = ^ < 5(j)n A 5vn > 

a=± n 

dz 

dp = , z^P± (3.45) 

z 

Note that the difference between this formula and the one in Section II. I is due to the 
choice t+ instead of t+ + t_ as the second variable. 

In the case of N x N matrix equations [dy — L,dt — A] = 0, we take the geometric 
data to consist of surfaces F with N punctures Pa, and E, Q to be real-normalized Abelian 
integrals with poles of order n and m at Pq. Let Oq, be given coefficients, with ai = 1. 
The leaf A4 corresponding to the above specifications for E, Q, combines with the space 
of parameters aa to a product space A4 x C^~^ of dimension g + N — 1. On the fibration 
j^g+N-i ai3oyg i]^Q product M. x C"^"^, we can construct the symplectic form 

g+N-l 

^M= Y. Q{ls)^E{^s) (3.46) 

s=l 

Now corresponding to the preceding geometric data are local coordinates Zq, near each 
puncture Pq, given in analogy with (3.21) and (3.22) by = a^^"", polynomials Qa^j^{z~^) 
which are the singular parts of Q near P^-, and thus a vector Baker-Akhiezer function 
'^{z\ X, y, t) = {'^a{z; X, y, t))^^i with the following essential singularity near Pp 

oo 

^a{z;x,y,t) = (5«/3 + X!^"^(^'2/'*)^^)^^P(^^'^ + «^^^"2/ + Q^,+(^^^)*) (3.47) 

s=l 

As in Section II. I, there exist then unique matrix operators L,A so that {dy — L)^ = 
{dt — ^)^' = 0. They have expressions of the form in (2.93) and subsequent equations. We 
observe that the case where Q has only simple poles (m^ = 1, Qp^+ = z^^) is the matrix 
generalization of the scalar case considered earlier, where Q still has the interpretation of 
a quasi-momentum. 
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On the space of such operators, we had defined in Theorem 8 a symplectic form u). 
Here again, the geometric and the formal symplectic forms (3.46) and (2.96) correspond 
to one another under the basic geometric correspondence Q. More precisely, 

N 

UJM = Y. ^^^Po. < **(^^^^ - L^^^SAra)-^ >dp = g^Um) (3.48) 
a=l 

where 

iUm = ReSocTr < %{A^^^SL - L^^^SA^) A S^^q > dk. (3.49) 
In the case = 1, we have A = dx, and this relation reduces to 

AT 

ujj^ = - ^ Resp^ < ^'*5L A 5* > dp = g^{uj) 

with CO given by (2.96). 

Similarly, our formalism can easily identify the action coordinates for the elliptic 
Calogero-Moser system, an issue which had been resolved only relatively recently [29] . We 
recall that the elliptic Calogero-Moser system is a system of N identical particles on a line, 
interacting with each other via the potential V{x) — p{x) 

Xi = 4j2p'{xi-Xj), Pi^) = ^^- (3-50) 

Here p{x) = p{x\lij,(jj') is the Weierstrass elliptic function with periods 2a;, 2a;', and a;, a;' 
are fixed parameters. The complete solution of the elliptic Calogero-Moser system was 
constructed by geometric methods in [35]. There an explicit Lax pair (L, M) was found, 
depending on a spectral parameter z varying on the torus C/(2wZ -|- 2a;'Z). Thus the 
dynamical system (3.50) is equivalent to the Lax equation L = [M,L], with L and M 
N X N matrices given by 

Lij{z) = piSij + 2(1 - Sij)^{xi - Xj,z), Pi = Xi 

Mij{z) = 25ij J2 Pi^i - ^j) + 2(1 - SijWixi - Xj) (3.51) 
and ^{x,z) is the function 

$(x,z) = 4^eC(^)^ (3.52) 

with a{z), C{z) the usual Weierstrass elliptic functions. In view of the Lax equation, the 
characteristic polynomial R{k,z) = det(2A; -|- L[z)) is time-independent, and defines a 
time-independent spectral curve F 

N 

R{k, z) = J2 n{z)k' = (3.53) 

i=0 
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where the ri{z) are eUiptic functions of z. The Jacobians of the spectral curves T are levels 
of the involutive integrals of the system. In particular, we obtain angle variables ipi by 
choosing 27r-periodic coordinates on them. However, as noted earlier, the identification of 
the canonically conjugate action variables is more difficult and has been carried out only 
recently [29]. In this work, it was shown that the Calogero-Moser sytem can be obtained 
through a Hamiltonian reduction of a Hitchin system, and as a result, the action- variables 
are the periods of the differential kdz along the A-cycles of the spectral curve V. 
We can derive this result directly from our approach. From our viewpoint, the leaf 
M. corresponding to the elliptic Calogero-Moser system is given by (F, k, z), where F is a 
Riemann surface of genus g — N, kisa, function with simple poles at g points Pi, ■ ■ ■ , Pg, dz 
is a holomorphic Abelian differential whose periods as well as integrals Jp" dz, 2 < a < g, 
are all on the lattice spanned by 2a; and 2a;', and the residues of kdz at Pa-, 2 < a < N 
are given by 

Resp„ (kdz) = 1, 2<a<N. (3.54) 

The analogues of the Baker-Akhiezer function and its dual are respectively the column 
vector C{P) = (ci, . . . , cjv) and the row vector C^[P) — (c]*", . . . , c^), satisfying 

{2k + L{z))C = C+{2k + L{z)) = 0, 

^ cM-^i, ■2) = 1, Yl "^tH^i, ^) = 1- (3-55) 

i i=l 

Here P — {k,z) is a point of the spectral curve F. The vectors C{P) and C^{P) are 
meromorphic functions on F outside the points Pq; on F corresponding to 2; = 0, and have 
each poles. We denote these poles by 71, • • • , and 7^, . . . , 7^^, respectively. Near the 
points Pa, these vectors have the form 

Ci{z) = {ct + 0{z))e^^'~\ c+iz) = {cf'+ + 0(2))e--^^" (3.56) 

where the coefficients cf satisfy 

c] = 1; ^c°= = for a > 1 

i 

c]'+ = 1; ""i'^ = for a > 1 (3.57) 

i 

The geometric symplectic form a;x constructed in Theorem 11 becomes in this case 



N , N 



e/V^ , / N , N 1 < A 5LC > , , , 

= S{2_^k{'ys)dz} = - 2^Resp„ dz, (3.57) 



2^ <C+C> 

s=l a=l 
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where < f'^g > denotes the usual pairing between column vectors and row vectors. Sub- 
stituting in the expansion (3.56), we obtain 



JV , TV 




(3.59) 



8—1 i 



This identifies our geometric symplectic form with the canonical symplectic form for the 
Calogero-Moser dynamical system. Since by construction (c.f.(3.42)), the geometric sym- 
plectic form admits the periods of kdz around Ai cycles as action variables dual to the 
angle variables on the Jacobian, our argument is complete. 



A. Non-linear WKB Methods in Soliton Theory 

We have seen that soliton equations exhibit a unique wealth of exact solutions. Nev- 
ertheless, it is desirable to enlarge the class of solutions further, to encompass broader 
data than just rapidly decreasing or quasi-periodic functions. Typical situations arising in 
practice can involve Heaviside-like boundary conditions in the space variable x, or slowly 
modulated waves which are not exact solutions, but can appear as such over a small scale 
in both space and time. 

The non- linear WKB method (or, as it is now also called, the Whitham method 
of averaging) is a generalization to the case of partial differential equations of the clas- 
sical Bogolyubov-Krylov method of averaging. This method is applicable to nonlinear 
equations which have a moduli space of exact solutions of the form uo{Ux + Wt + Z\I). 
Here uo{zi, . . . , Zg\I) is a periodic function of the variables zf, U = (t/i, . . . , Ug), W = 
{Wi, . . . , Wg) are vectors which like u itself, depend on the parameters / = (/i, . . . , In), 
i.e. U — U{I), V — V{I). (A helpful example is provided by the solutions (3.10) of the 
KP equation, where / is the moduli of a Riemann surface, and U, V, W are the S^-periods 
of its normalized differentials dOi, d02, and dOa.) These exact solutions can be used as a 
leading term for the construction of asymptotic solutions 



u{x, t) = UQ{e-^S{X, T) + Z{X, T)\I{X, T)) + eui{x, t) + e^U2{x, t) + • • • , (4.1) 



where / depend on the slow variables X = ex,T = et and and e is a small parameter. If 
the vector-valued function S{X, T) is defined by the equations 



then the leading term of (4.1) satisfies the original equation up to first order one in e. 
All the other terms of the asymptotic series (4.1) are obtained from the non-homogeneous 



IV. WHITHAM EQUATIONS 



dxS = C/(/(X, T)) = U{X, T), drS = W{I{X, T)) = W{X, T), 



(4.2) 
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linear equations with a homogeneous part which is just the hnearization of the original 
non-linear equation on the background of the exact solution uq. In general, the asymptotic 
series becomes unreliable on scales of the original variables x and t of order e~^. In order to 
have a reliable approximation, one needs to require a special dependence of the parameters 
I{X, T). Geometrically, we note that e~^S{X, T) agrees to first order with Ux + Vy + Wt, 
and X, y, t are the fast variables. Thus u{x, t) describes a motion which is to first order the 
original fast periodic motion on the Jacobian, combined with a slow drift on the moduli 
space of exact solutions. The equations which describe this drift are in general called 
Whitham equations, although there is no systematic scheme to obtain them. 

One approach for obtaining these equations in the case when the original equation is 
Hamiltonian is to consider the Whitham equations as also Hamiltonian, with the Hamil- 
tonian function being defined by the average of the original one. In the case when the 
phase dimension g is bigger than one, this approach does not provide a complete set of 
equations. If the original equation has a number of integrals one may try to get the com- 
plete set of equations by averaging all of them. This approach was used in [62] where 
Whitham equations were postulated for the finite-gap solutions of the KdV equation. The 
geometric meaning of these equations was clarified in [26]. The Hamiltonian approach 
for the Whitham equations of (1-|-1) -dimensional systems was developed in [23] where the 
corresponding bibliography can also be found. 

In [36] a general approach for the construction of Whitham equations for finite-gap 
solutions of soliton equations was proposed. It is instructive enough to present it in case 
of the zero curvature equation (2.1) with scalar operators. 

Recall from Sections III. A and III.B that the coefficients Ui{x,y,t), Vj{x,y,t) of the 
finite-gap operators Lq and Aq satisfying (2.1) are of the form (c.f. (3.10)) 

Ui^Ui^o{Ux + Vt + Wt + Z\I), Vj ^Vj^o{Ux + Vt + Wt + Z\I), (4.3) 

where and Vj^o are differential polynomials in ^-functions and / is any coordinate 
system on the moduli space A4 g{n, m). 

We would like to construct operator solutions of (2.1) of the form 

L = Lo + eLi + ■ ■ ■ , A = Ao + eAi + • • • , (4.4) 

where the coefficients of the leading terms have the form 

ui = w,,o(e-^5(x, y, T) + z{x, y, T) |/(x, y, t)), 

= vj,o{e-'S{X, Y, T) + Z(X, F, T) |/(X, F, T)) (4.5) 

From Section III.B, we also recall that J\fg{n,m) is the bundle over M.g{n,m) with the 
corresponding curve F as fiber. If / is a system of coordinates on M.g{n, m), then we may 
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introduce a system of coordinates {z, I) on J\fg{n,m) by choosing a coordinate along the 
fiber r. The Abehan integrals p, E, Q are multi- valued functions of (A, /), i.e. p = p{X, /), 
E = E{X,I), Q = Q{X,I). If we describe a drift on the moduli space of exact solutions by 
a map (X, F, T) — > / = I{X, Y, T), then the Abelian integrals p, E, Q become functions of 
z, X, Y, T, via 

p{z,X,Y,T)=p{z,I{X,Y,T)), 
E{z,X,Y,T) = E{z,I{X,Y,T)), 
Q{z,X,Y,T) = Q{z,I{X,Y,T)). 

The following was established in [36] : 

Theorem 12. A necessary condition for the existence of the asymptotic solution (4-4) 
with leading term (4-3) and bounded terms Li and Ai is that the equation 

dpfdE dQ\ dEfdp dQ\ dQ { dp dE\ ^ 

+ ^ ^ - ^ = (4-6) 



dz \dT dY ) dz \dT dX ) dz \ dY dX 

be satisfied. 

The equation (4.6) is called the Whitham equation for (2.1). It can be viewed as a 
generalized dynamical system on Aig{n,m), i.e., a map {X,Y,T) — > M.g{n,m). Some of 
its important features are: 

• Even though the original two-dimensional system may depend on y, Whitham so- 
lutions which are y-independent are still useful. As we shall see later, this particular case 
has deep connections with topological field theories. If we choose the local coordinate z 
along the fiber as z = E, then the equation simplifies in this case to 



OtP = dxQ (4.7) 
We note that it followed immediately from the consistency of (4.2) that we must have 




Thus (4.7) is a strengthening of this requirement which encodes also the solvability term 
by term of the linearized equations for all the successive terms of the asymptotic series 
(4.3). 

• Naively, the Whitham equation seems to impose an infinite set of conditions, since 
it is required to hold at every point of the fiber F. However, the functions involved are 
all Abelian integrals, and their equality over the whole of F can actually be reduced to a 
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finite set of conditions. To illustrate this point, we consider the Y-independent Whitham 
equation on the moduli space of curves of the form 



= 1[{E - Ei) = R{E) 



i=l 

Then the differentials dp and dQ are given by 

rEa EdE 
dp=^^{E-^^^) 

We view p and Q as functions of {E; Ei), with E the coordinate in the fiber F, and the 
coordinates on the moduli space of curves. Near each branch point E'j, y/E — Ei is a local 
coordinate and we may expand 

p = p{Ei) + ay/E - Ei + 0{E - E^) 

Q = Q{Ei) + (3^JE -Ei + 0{E - Ei) (4.8) 

Differentiating with respect to X and T, keeping E fixed, we find that the leading singular- 
ities of drp and dxQ are respectively -^--^^=dTEi and -^-j^==dxEi. Since f = 
we see that the equation (4.7) implies 

dTE,= {^)\^dxE, (4.9) 

Conversely, if the equation (4.9) is satisfied, then drp—dxQ is regular and normalized, and 
hence must vanish. Thus the equation (4.7) is actually equivalent to the set of differential 
equations (4.9). 

• The equation (4.7) can be represented in a manifestly invariant form, without explicit 
reference to any local coordinate system z. Given a map {X,Y,T) — > A4 g{n,m), the pull 
back of the bundle jVg (n, m) defines a bundle over a space with coordinates X, Y, T. The 
total space jV^ of the last bundle is 4-dimensional. Let us introduce on it the one-form 

a = pdX + EdY + QdT, (4.10) 

Then (4.7) is equivalent to the condition that the wedge product of da with itself be zero 
(as a 4-form on M.'^) 

da A da = 0. (4.11) 
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• It is instructive to present the Whitham equation (4.7) in yet another form. Because 
(4.7) is invariant with respect to a change of focal coordinate we may use p = p{z, I) by 
itself as a local coordinate. With this choice we may view E and Q as functions of X, Y 
and T, i.e. E = E{p,X,Y,T), Q = Q{p,X,Y,T). With this choice of local coordinate 
(4.7) takes the form 

dTE-dYQ + {E,Q} = 0, (4. 12) 

where {•, •} stands for the usual Poisson bracket of two functions of the variables p and X, 
i.e. 

if, 9} = fp9x-9pfx- 

• In Theorem 12 we had focused on constructing an asymptotic solution for a single 
equation. This corresponds to a choice of A, and thus of an Abelian differential Q, and 
the Whitham equation is an equation for maps from {X,Y,T) to Aig{n,m). As in the 
case of the KP and other hierarchies, we can also consider a whole hierarchy of Whitham 
equations. This means that the Abelian integral Q is replaced by the really normalized 
Abelian integral Qi which has the following form 

Qi = K' + 0{K-^), K'' = E, (4.13) 

in a neighborhood of the puncture P (compare with (2.48)). The result is a hierarchy of 
equations on maps of the form 

(X = Ti, F = T2, T = T3, T4, • • •) ^ Mg{n). (4.13) 

The whole hierarchy may be written in the form (4.11) where we set now 

a = J2^i'iTi. (4.14) 

i 

B. Exact Solutions of Whitham Equations 

In [38] a construction of exact solutions to the Whitham equations (4.7) was proposed. 
In the following section, we shall present the most important special case of this construc- 
tion, which is also of interest to topological field theories and supersymmetric gauge the- 
ories. It should be emphasized that for these applications, the definition of the hierarchy 
should be slightly changed. Namely, the Whitham equations describing modulated waves 
in soliton theory are equations for Abelian differentials with a real normalization (3.29). 
In what follows we shall consider the same equations, but where the real-normalized differ- 
entials are replaced by differentials with the complex normalization (3.30). As discussed 
in Section III. A, the two types of normalization coincide on the subspace corresponding 
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to M-curves, which is essentially the space where all solutions are regular and where the 
averaging procedure is easily implemented. Thus the two forms of the Whitham hierarchy 
can be considered as different extensions of the same hierarchy. The second one is an 
analytic theory, and we shall henceforth concentrate on it. 

In this subsection and in the rest of the paper, we shall restrict ourselves to the 
hierarchy of "algebraic geometric solutions" of Whitham equations, that is, solutions of 
the following stronger version of the equations (4.11) 

^E = {n„E}, (4.15) 

We note that the original Whitham equations can actually be interpreted as consistency 
conditions for the existence of an E satisfying (4.15). Furthermore, the solutions of (4.15) 
can be viewed in a sense as "y-independent" solutions of Whitham equations, since the 
equation (4.12) reduces to drE + {E,Q} = for F-independent solutions. They play 
the same role as Lax equations in the theory of (2+l)-dimensional soliton equations. As 
stressed earlier, y-indepcndcnt solutions of the Whitham hierarchy can be considered even 
for two-dimensional systems where the y-dependence is non-trivial in general. 

Our first step is to show that (4.15) defines a system of commuting flows on M.g{n). 
For the sake of simplicity, we assume that there is only one puncture. Let us start with 
a more detailed description of this space which is a complex manifold with only orbifold 
singularities. Its complex dimension is equal to (c.f. (3.31)) dim.Aig{n) = 4g + n — 1, 
and we had constructed a set of local coordinates for it in Theorem 10 and subsequent 
discussion. Here we require the following slightly different set of coordinates (details can 
be found in [38]). The first 2g coordinates are still the periods of d£^, 

TAuE = / dE, TB„E = (f dE (4.16) 

The differential dE has 2g -\- n — 1 zeros (counting multiplicities). When all zeroes are 
simple, we can supplement (4.16) by the 2g -\- n — 1 critical values Eg of the Abelian 
differential E, i.e. 

E, = E{q,), dE{q,)^0, s = 1, . . . , 2^7 + n - 1, (4.17) 

In general, dE may have multiple zeroes, and we let D = X] A*s9s be the zero divisor of dE. 
The degree of this divisor is equal to J^^^ //g = 2^ — 1 -|- n. Consider a small neighborhood 
of Qs, viewed as a point of the fibration Afg{n), above the original data point mo in the 
moduli space M.g{n). Viewed as a function on the fibration, E is a deformation of its 
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value E{z,mo) above the original data point, with multiple critical points Qg. Therefore, 
on each of the corresponding curve, there exists a local coordinate Wg such that 

E = w^^+\z,m)+ J2 Es,,{m)wl{z,m). (4.18) 

The coefficients Es^i{m) of the polynomial (4.18) are well-defined functions on Mg{n). 
Together with T^j,s, 'TBi,E, they define a system of local coordinates on A^g(n).If //g = 1, 
then Eg^Q clearly coincides with the critical value E{qg) from (4.17). 

Let V be the open set in M.g{n) where the zero divisors of dE and dp, namely the 
sets {z\dE{z) — 0} and {z]dp{z) — 0}, do not intersect and let be the open set in 
M.g{n) where all zeros of dE are simple. 

Theorem 13. The Whitham equations (4-15) define a system of commuting meromorphic 
vector fields (flows) on Aig{n) which are holomorphic on V C M.g{n). On the open set 
V' nV, the equations (4-15) have the form 

d d 

rA,,E^O, TB^^E^O, (4.18) 



dTj " ' dTj 



dT,Eg = ^{qg)dxEg. (4.19) 



Indeed, the equations (4.15) are fulfilled at each point of F, and thus 

A£ds^ ^^wa. (£dE) - f (£d«.) (4.20) 

The functions dE/dp and dfli/dp are linearly independent. It follows that the periods of 
dE arc constants. The equations (4.19) coincide with (4.15) at the point qg (where dE 
equals zero). 

In order to complete the proof of Theorem 13, it suffices to show that (4.18-4.19) 
imply (4.15). The equation (4.19) implies that the difference between the left and right 
hand sides of (4.15) is a meromorphic function f{z) on F. This function is holomorphic 
outside the puncture P and the zeros of dp. At the puncture P, the function f{z) has a 
pole of order less or equal to (n — 2). However, f{z) equals zero at zeros of dE. Hence, the 
function g{z) = fiz)^ is holomorphic on F and equals zero at P. Therefore, f{z) = 
identically. Theorem 13 is proved. 

An important consequence of Theorem 13 is that the space M.g{n) admits a natural 
foliation, namely by the joint level sets of the functions TAi,E, TBi,E, which are smooth 
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{2g + n — 1) -dimensional submanifolds, and which are invariant under the flows of the 
Whitham hierarchy (4.15). We shaU sometimes refer to the leaves M. of this foliation as 
large leaves^ to stress their distinction from the ^-dimensional leaves Ai of the canonical 
foliation of M.g{n,m). 

The special case of the construction of exact solutions to (4.15) in [38] may now be 
described as follows: the moduli space M.g{n, m) provides the solutions of the first n + m- 
flows of (4.15) parametrized by 3g constants, which are the set of the last three coordinates 
(3.38-3.39) on Mg{n,m). 

Theorem 14. Let Ti, i = l,...,n + m, TAi,E, rBi,E,rAi,Q, TAi,Q, ai be the coordinates 
on A4g{n,m) defined in Theorem 10. Then the projection 

Mg{n, m) Mg{n) 
{T,E,Q)^{T,E) (4.21) 

defines (F, E) as a function of the coordinates on A4g{n, m). For each fixed set of param- 
eters TAi,E, rBi,E,rAi,Q, TAi,Q, ai, the map {TiY^I'^'^ Mg{n) satisfies the Whitham 
equations (4-- 15). 

Proof. Let us use E[z) as a local coordinate on F. Then as we saw earlier, the equations 
(4.15) are equivalent to the equations dTiP{E,T) = dx^i{E,T). These are the compati- 
bility conditions for the existence of a generating function for all the Abelian differentials 
dfli. In fact, if we set 

dA = QdE (4.22) 
then it follows from the definition of the coordinates that 

dT,dX = dili, dxdX = dp, (4.23) 

(For the proof of (4.23), it is enough to check that the right and the left hand sides of it 
have the same analytical properties.) 

Theorem 15. We consider the same parametrization of A4g{n,m) as in Theorem 14- 
Then as a function of the parameters Ti, 1 < i < n + m, the second Abelian integral 
Q(j),T) satisfies the same equations as E, i.e. 

dT,Q = {ni,Q}. (4.24) 

Furthermore 

{E,Q} = 1. (4.25) 

We note that (4.25) can be viewed as a Whitham version of the so-called string 
equation (or Virasoro constraints) in a non-perturbative theory of 2-d gravity [19] [66]. 
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C. The T-Function of the Whitham Hierarchy 

The solution of the Whitham hierarchy can be succinctly summarized in a single 
T-function defined as follows. The key underlying idea is that suitable submanifolds of 
A4g{n, m) can be parametrized hy 2g -\- N — 1 -\- X]^=i(^a + "^a) Whitham times T^, to 
each of which is associated a "dual" time Tda, and an Abelian differential dVtA- We begin 
by discussing the simpler case of one puncture, = 1. Recall that the coefficients of the 
pole of dA has provided n + m Whitham times Tj — — j'Rjes{z^ dX) . Their "dual variables" 
are 

Tdj =Res{z-^dX), (4.26) 

and the associated Abelian differential are the familiar dQi of (3.4) (complex normalized). 
When g > 0, the moduli space Aig{n, m) has in addition 5g more parameters. We consider 
only the foliations for which the following 3g parameters are fixed 

/ dE, (f dE, (f dQ. (4.27) 

JAk JBu JAk 

Thus the case ^ > leads to two more sets of g Whitham times each 

ak= i dA, = (f dQ (4.28) 
Jau Jbu 



Their dual variables are 

ClDk 



2m 



j dA, TEk = ^.j EdX (4.29) 



The corresponding Abelian differentials are respectively the holomorphic differentials du^ 
and the differentials dO^", defined to be holomorphic everywhere on F except along the Aj 
cycles, where they have discontinuities 

dl]f + - dJlf- = 5jkdE (4.30) 

We denote the collection of all 2g + n + m times by Ta = {Tj^ak.T^). In the case of 
A^ > 1 punctures, we have 2g + X]a(^" + TTia) times (Tq,j, a/j, T^) and 3A^ — 3 additional 
parameters for A4g{n,m), namely the residues of dQ, dE, and dA at the punctures Pq, 
2 < q; < A^ (c.f. (3.35-3.36)). For simplicity, we only consider the leaves of Mg{n,m) 
where 

Resp„(dQ) = 0, Resp„(d£;) = fixed, 2<a<N, (4.31) 

and incorporate among the Ta the residues of dA at Pq, 2 < ct < A^, = Resp^(dA) 
(c.f. (3.35)). The dual parameters to these A^ — 1 additional Whitham times are then the 
regularized integrals 

Roa = - / dA, 2 < a < AT. (4.32) 
Jpi 
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More precisely, recall that the Abelian integral A has been fixed by the condition that its 
expansion near Pi, in terms of the local coordinate Zq, defined by E, have no constant 
term. Near each Pq,, 2 < ck < AT, in view of (4.31), it admits an expansion of the form 



Tla.~\~'^a. rji 



K^a)= ^+Resp„(dA)log2;a + A« + 0(2;a) (4.33) 

For each a, 2 < a < A^, we define the right hand side of (4.32) to be Aq. The Abelian 
differential dOa^"* associated with R'^ is the Abelian differential of the third kind, with 
simple poles at Pi and Pa-, and residue 1 at P^. In summary, we have the following table 

Times Dual Times Differentials 



iRespJ^^dA) RespJ^-MA) dfi^ 

f^.dA -^f^.dA do;; (4-34) 



■ dA dwfc 
Resp^dA -/^i"dA doi^^ 

§B, C dnf 



We can now define the r-function of the Whitham hierarchy by 

t(T) = e^(^) 

1 1 ^ 

•^(^) = o E ^^^DA + r- E (^kT^^iAk n Bk) (4.35) 

A k=l 

where AkOBk is the point of intersection of the and Bf, cycles. In the case Resp^d£? = 
we have then (see [38]) 

Theorem 16. The derivatives of T with respect to the 2g+^^{na+ma) + N—l Whitham 
times Ta are given by 



1 ^ 

Qt.T ^Tda + — Y1 ^a„ATi^E{Ak n Bk) 
fc=i 

^T„,,,T^,,-^ = Resp„ (4d^^/3,j) 
^(-.^)-^=2^£ ^^^^ 



-77 /'dOAdnBdnc\ . . 

sisc^ = E ^^^^'^^ ( aEdQ I ^^-^^^ 
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These formulae require some modifications when Resp^dE^ ^ 0, which is a case of 
particular interest in supersymmetric QCD (see Section VI below). In particular, the first 
derivatives with respect to become [17] 

1 - 

where c^^f^ is an integer which is antisymmetric in a and /9. It is then easy to see that the 
second derivatives are modified accordingly by constant shifts, while the third derivatives 
remain unchanged. 

We would like to point out that the 1g + ^J^a. + "^a) -\- N — \ submanifolds of 
M.g{n^m) defined by fixing (4.27) as well as the residues of di? and dQ are yet another 
version of the 2^ + n — 1 large leaf M. of the foliation of Jv[g{n) encountered earlier in 
the case of one puncture. Indeed, imbedding JV[g{n) in Aig{n, 1) by choosing = k 
would fix two Whitham times and T^+i, as we saw after Theorem 10. This reduces the 
dimension 2g + n + 1 to the desired dimension 2g + n — 1. 

We observe that the first derivatives of J- give the coefficients of the Laurent expan- 
sions and the periods of the form A. The second derivatives give the coefficients of the 
expansions and the periods of the differentials dfi^. In that sense the function encodes 
all the information on the Whitham hierarchy. The formulae for the first and the second 
derivatives with respect to the variables are the analogues in the case of the averaged 
equations of the corresponding formulae for the r-function of the KP hierarchy. The formu- 
lae for the third derivatives are specific to the Whitham theory. As we shall see later, they 
are reminiscent rather of marginal deformations of topological or conformal field theories 
and of special geometry. 

Finally, it may be worth noting that the expression (4.35) for can be elegantly 
summarized as 

T = J dAAdA. (4.37) 

V. TOPOLOGICAL LANDAU-GINZBURG MODELS 
ON A RIEMANN SURFACE 

In this section, we shall show that each 2g + n — 1 leaf of the foliation of Mg{n) (or, 
equivalently, of the foliation of Aig{n, 1) upon imbedding Jlig{n) into J^g{n, 1)) actually 
parametrizes the marginal deformations of a topological field theory on a surface of genus 
g. Furthermore, the free energy of these theories coincides with the restriction to the leaf 
of the exponential of the r-function for the Whitham hierarchy. We begin with a brief 
discussion of some key features of topological field theories [63-65] [11-13]. 
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A. Topological Field Theories 

In general, a two-dimensional quantum field theory is specified by the correlation 
functions < 0(^i) ■ ■ ■(/>(2;jv) >g of its local physical observables (j)i{z) on any surface F of 
genus g. Here 4>i{z) are operator- valued tensors on V. The operators act on a Hilbert 
space of states with a designated vacuum state \Vt >. The correlators < (t>{zi) ■ ■ -(^{zn) > 
usually depend on the background metric on F and on the location of the insertion points 
Zi- In particular, they may develop singularities as Zi approaches Zj. Equivalently, the 
operator product (f)i{zi)4)j{zj) may develop singularities. For example, in a conformal field 
theory, (f)i{zi)(f)j{zj) will have an operator product expansion of the form 

4>i{zi)4>j{zj) = - Zj)^^^^^~^'"4>k{zj) + descendants (5.1) 

k 

where hi is the conformal dimension of If we let (I)q{z) be the field corresponding to 
|f2 >, under the usual states <-> fields correspondence, then we obtain a metric by setting 

= 4 (5.2) 

Using 'qij to raise and lower indices, and noting that < 4>o[z) >o= 1, we can easily recognize 
Cijk — CijTjki as the three-point function on the sphere 

Cijk =< (l>iizi)(j)j{zj)(t)k{zk) >o, (5.3) 

which is actually independent of the insertion points Zi, Zj, Zk by SL(2,C) invariance. 

Topological field theories are theories where the correlation functions are actually 
independent of the insertion points Zi. Thus they depend only on the labels of the fields 
and the genus g of F. This independence implies that the OPE of (5.1) contains no 
singularity, and that for all practical purposes, the operator product (f)i{zi)(l)j{zj) can be 
replaced by the formal operator algebra 

k 

The associativity of operator compositions translates into the associativity of the operator 
algebra (5.4). Furthermore, the operator algebra is commutative, since factorization of 
the 4-point function in the s and t channels must give the same answer. If we assume 
that in a topological field theory, physical correlation functions can be factored through 
only physical states (as is the case, if the Hilbert space of physical states arises as the 
cohomology of a nilpotent BRST operator Q acting on a larger Hilbert space containing 
spurious states), then the correlation functions of the theory on surfaces of all genera can 
all be expressed explicitly, through factorization, in terms of the structure constants Cijk- 
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B. Deformations of a Topological Field Theory 

We shall be particularly interested in the the case where the topological field theory 
arises as a deformation by parameters ti of a fixed theory. In a topological field the- 
ory, the physical fields (pi have clearly scaling dimension 0. However, in presence of a 
BRST symmetry, they actually fit in a multiplet $i = {(pf'^ — (pi^ (f)\^^), where the 
"descendants" are tensors of scaling dimension d which satisfy the descent equations 
ci(^(^) = {(5, (/>*-'^''~^-*}. In particular, (^p"* can be integrated on the surface F. If we let /(O) 
be the Lagrangian of the original theory, we can deform the theory to another theory with 
Lagrangian 

m = m-J2t^ 1 4>?^ (5.5) 

The structure constants Cijk{t) of the new theory are then given by 

Cijk{t) =< (j)i(j)j(j)kex.p{y2u / (Pi^^) >o (5.6) 

They define in turn a topological field theory. For deformations around a topological 
conformal field theory, i.e., a topological field theory whose stress tensor is already traceless 
before restricting to physical states, the structure constants Cijk{t) are known to satisfy 
the key compatibility condition diCijk = dkCiji [11-13]. This means that we can find a 
function J^, called the free energy and formally denoted by 

T{t) =< exp(^ti f ^i) >o, (5.7) 

i 

which satisfies the third derivative condition 

dfjkm = Cijk (5.8) 

In terms of the commutativity and associativity conditions of the operator algebra with 
structure constants Cijk{t) become a system of differential equations of third order, called 
the Witten-Dijkgraaf-Verlinde-Verlinde (WDVV) equation. 

C. The Framework of Solitons 

In their original work [11-12], Dijkgraaf, Verlinde, and Verlinde derived an explicit 
expression for the free energy jF(t) in the case of topological Landau-Ginzburg theories. 
These are the topological theories arising from twisting an N=2 superconformal field theory, 
which is itself obtained by following the renormalization group flow to the fixed point of 
a Landau-Ginzburg model. Although the renormalization group flow modifies the kinetic 
terms in the Landau-Ginzburg action, the superpotential W{x) remains unchanged, and 
thus characterizes both the associated superconformal and topological models [30] [61] [67]. 
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Our goal in this section is to show how this theory can fit in the framework of sohtons, 
and to exhibit the natural emergence of the difi^erential QdE and Whitham times. 

We consider first the case of genus 0, with F = {z G CUoo}. The role of the 
superpotential W{x) is played in our context by the Abelian integral E with a unique pole 
of order n at oo. We consider then a leaf in the space A4o{n,l), characterized by the 
condition that Q = z, and E is of the form 

n-2 

E ^ z"" + +0{z-^) (5.9) 

This is of dimension n — 1, and can be parametrized by the n — 1 Whitham times Ta, 
A = 1, ■ ■ ■ , n — 1, with the other times fixed to = 0, T„+i = At each point E on 

the leaf, the primary fields (pi can be identified with dQi/dQ. The structure constants are 
defined by 

dncjidriciK , dnc,-)drici-\dficii.v 

< <P,<Pj >= Res(— ^1^), < <P,<Pj<Pk >= Res( ^ ^ (5-10) 

We note that dQ corresponds to the field defined by the vacuum. The following can be 
derived from Theorem 16: 

Theorem 17. LetJ^iT) he the r -function of the (complex normalized) Whitham hierarchy, 
restricted to the {n — 1) dimensional leaf described above. Then 
(i) the fields (pi anti-diagonalize the pairing rjij, i.e., < (pi(pj >= 5iJ^j^n', 
(a) J-'{T) is the free energy of the theory, i.e., f^T.T,Tfe"^(-^) ~ particular, T{T) 

satisfies the system of WDVV equations. 

More generally, the case of V of genus g (with one puncture to simplify the notation) 
has been treated in [22] [38]. In this case, the relevant leaf within A4f,(n, 1) is of dimension 
n — 1 + 2g and is given by the constraints 

Tn = 0, Tn+i 



n+1 



i dE = 0, f dE = fixed, 

JAk JBk 

/ dQ = (5.11) 

Thus the leaf is parametrized by the (n — 1) Whitham times T^, A = l,---,n — 1, and 
by the periods aj and of (4.28). The fields (pA of the theory need to be augmented 
accordingly. We take the 2g additional ones to be given by dujj/dQ and dfi^/dQ, where 
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the differentials da;^ and dHf are the ones associated to aj and , as described ear her 
in (4.30). 

Theorem 18. Let tjab and cabc be defined as 

Qs 

E dn^dOfidnc 
^'^^ dEdQ ' ^^-^^^ 

TO^/i gs zeroes ofdE, and the indices A,B,C running this time through the augmented 
set of n — 1 + 2g indices given by Ta = iTi, aj, Tej). Then 

(i) r]ij = Si+j^n, r]a^,{E,k) = 5j,k- All other pairings vanish; 

(ii) Let T{Ti,aj^TE,j) be the r-function of the Whitham hierarchy restricted to the leaf 
(5.11). Then d'^BC-^i'^) ~ cabc, where A runs through the augmented set of n — 1 + 2g 
indices. 

Note that in the genus case when Q — z, the sum in (5.12-5.13) over the residues at 
the zeros of dE reduces to the residue at infinity. 

Remarkably, the larger spaces M.g{n,m) can accommodate the gravitational descen- 
dants of the fields 0^. More precisely, consider for ^ = the leaf of the space M.o{n, mn+1) 
given by the following evident modification of the earlier normalization (5.11) 

Tin = 0, i= l,...,m, Tnm+i = • (5-14) 

nm + 1 

The space of Whitham times is automatically increased to the correct number by taking 
all the coefficients of QdE. The additional m{n — 1) fields may be identified with the first 
m gravitational descendants of the primary fields. Namely, the p-th descendant o'p{(f)i) of 
the primary field 0i is just This statement is a direct corollary of the following 

theorem. 

Theorem 19. The correlation functions given by < 4>a4>b4>c >— dABC-^ ^p{(pi) = 
dQi-^-pn/dQ satisfy the factorization properties for descendant fields 

< (Tp{^i)^B^C > = < crp-ii4>i)(l>j > < 4>k4>B4>c > (5.15) 

where 4)i,i = 1, . . . ,n — 1 are primary fields, and 4) a are all fields (including descendants). 

Factorization properties for descendant fields were derived by Witten [63-66]. We 
note that the completeness of the operator algebra requires a larger set of fields that just 
g, which is the dimension of the small leaves of the canonical foliation of M.g{n, m), and 



63 



which will be shown in the next section to be the dimension of the moduli space of vacua of 
certain supersymmetric gauge theories. This is one of the difficulties in establishing direct 
contact between topological field theories and supersymmetric gauge theories, although 
there has been progress in this direction [6] [44] [45] . 

VI. SEIBERG-WITTEN SOLUTIONS OF N=2 SUPERSYMMETRIC 

GAUGE THEORIES 

Moduli spaces of geometric structures are appearing increasingly frequently as the key 
to the physics of certain supersymmetric gauge or string theories. One recurrent feature 
is a moduli space of degenerate vacua in the physical theory. The physics of the theory is 
then encoded in a Kahler geometry on the space of vacua, or, in presence of powerful con- 
straints such as N=2 supersymmetry, in an even more restrictive special geometry, where 
the Kahler potential is dictated by single holomorphic function J-', called the prepotential. 
This was the case for Type IIA and Type IIB strings, where the vacua corresponding to 
compactifications on Calabi-Yau threefolds [31] [68] produce effective N=2 four dimensional 
supergravity theories. The massless scalars of such theories (in this case, the moduli of 
the Calabi-Yau threefold) must parametrize a manifold equipped with special geometry 
[9] [51] [57]. More recently, a similar phenomenon has been brought to light by Seiberg- 
Witten [52] [53] for N=2 supersymmetric gauge theories. Remarkably, the space of vacua 
of these theories, which is classically just a space of diagonalizable and traceless matrices, 
becomes upon quantization a moduli space of Riemann surfaces. The prepotential for 
the quantum effective theory can then be derived from a meromorphic one-form dA on 
each surface. A particularly striking feature of these effective theories, noticed by many 
authors [24] [28] [42-43] [49] [56] , is a strong but as yet ill- understood similarity with the 
Whitham theory of solitons. Indeed, the quantum spaces of vacua for many N=2 SUSY 
theories actually coincide with certain leaves of the canonical foliation of M.g{n,m) [39], 
the Seiberg-Witten form dA with the one-form QdE central to Whitham theory, and the 
effective prepotential of the gauge theories with the exponential of the r-function of the 
Whitham theory! The purpose of this section is to review some of these developments. 

A. N=2 Supersymmetric Gauge Theories 

We begin with a brief account of N=2 SUSY Yang- Mills theories in four dimensions 
with gauge group G [2]. The Yang-Mills gauge field A = Afj^dx^ is imbedded in an N=2 
gauge multiplet consisting of A, left and right Weyl spinors Xl and Xr, and a complex 
scalar field (p, with all fields valued in the adjoint representation of G. The requirement of 
N=2 SUSY and renormalizability fixes uniquely the action 
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where g is the couphng constant, 6 is the instanton angle, and we have written exphcitly 
only the bosonic part of the action. The classical vacua are given by the critical points 
of the action. In this case, they work out to be ^ = 0, is constant (up to a gauge 
transformation), and 

[<^,0t]=o (6.2) 

Thus (j) must lie in the Cartan subalgebra. For G=SU(A^c) (-^c is commonly referred to as 
number of "colors"), we set 



/ai \ 







02 



J]afe = 0. (6.3) 



fc=i 



7 



Thus the classical space of vacua is parametrized by the Ofc, up to a Weyl permutation. 

For a generic configuration a^, we have aj ^ for any j ^ A;, and the gauge group 
SU(A^c) is spontaneously broken down to U(l)'^'' ^ . At the quantum level, we expect then 
the space of inequivalent vacua to be parametrized by N^. — 1 parameters ajt (thought of as 
renormalizations of the a^, = 0), with each vacuum corresponding to a theory of 

Nc — 1 interacting U(l) gauge fields Aj, i.e., Nc — l copies of electromagnetism. In the weak 
coupling regime, we expect singularities at ak = ai, where the gauge symmetry is suddenly 
enhanced. Since N = 2 SUSY remains unbroken, each gauge field Aj is part of an N=2 
SUSY U(l) gauge multiplet {Aj,XLj.,XRj.,(f)j)., all in the adjoint representation of U(l). 
Again, the action for a theory with such a field content is fixed by N=2 supersymmetry. 
To leading order in the low momentum expansion for these fields, it must be of the form 

le// = — / d^a; r(ImT^'=)F. A *Fk + {ReT^'')Fj A + d^'^^ A d^f^l + fermions (6.4) 

Stt Jm4. 

where 

for a suitable complex and analytic function jF(a, A), called the prepotential. We note that 
the prepotential .7-" is a function not just of the vacua parameters a^, but also of a scale A 
introduced by renormalization. 

Thus the physics of the quantum theory is encoded in a single function J^. What is 
known about J-"! To insure the positivity of the kinetic energy, we must have 

Im— — — > 0. (6.6) 
oajOak 
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Geometrically, T defines then a Kahler metric on the quantum moduli space by 



Furthermore, at weak-coupling A << 1, ^ can be evaluated in perturbation theory. For 
pure SU(Arc) Yang-Mills, one finds 

^(«.^)=|^i:'4-8i. EK-^.O^'^g^^^^i^^ + E^^A-." (6.8) 

fe=l j,fc=l d=l 

The first term on the right hand side is the classical prepotential. The second term is the 
perturbative one-loop quantum correction. In view of N=2 non-renormalization theorems, 
it is known that higher loops do not contribute. The third term is the instanton contri- 
bution, consisting of d-instanton processes for all d. We observe that the expansion (6.8) 
implies in particular that T has non-trivial monodromy around = in the A << 1 
regime. The exact solution of N=2 Yang-Mills theories is reduced in this way to finding a 
holomorphic T satisfying the constraints (6.6) and (6.8). 

We have just described the main problem for N=2 SUSY pure SU(A^) Yang-Mills 
theories. However, the same problem should be addressed for general N=2 SUSY gauge 
theories with gauge group G, with matter fields ( "hypermultiplets" ) in a representation R 
of G. As in the case of pure Yang-Mills, the Wilson effective Lagrangian of these theories is 
dictated by a prepotential J-'cRia, A), and the problem is to determine J-'g,r{0', A), subject 
to the constraints (6.6) and (6.8), where the right hand side of (6.8) has been modified 
to incorporate the contributions of the hypermultiplets. For example, in presence of Nf 
hypermultiplets in the fundamental representation of bare masses rrii, 1 < i < Nf, the 
one-loop correction to the prepotential for the SU(A^c) theory contains the additional term 

2^ 2^(afc + m,) log J2^- 

B. The Seiberg-Witten Ansatz 

The requirements that T have monodromy and a Hessian with positive definite imag- 
inary part, suggest an underlying non-trivial geometry on the quantum space of vacua. In 
[52] [53] , Seiberg and Witten made the fundamental Ansatz that 

• For each A, the quantum moduli space should parametrize a family of Riemann sur- 
faces r(a. A) of genus g = — 1, now known as the spectral curves of the theory; 

• on each r(a,A), there is a meromorphic one-form dA; 
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• ^ is determined by the periods of dA * 





(6.9) 



= 0'D,k 



The gauge theories under consideration contain dyons, i.e., particles which carry both 
electric and magnetic charges. Let (n, m) e Z-'^^"^ x Ti^''~^ be their charges, with (n^, m^) 
the charge with respect to the z-th U(l) factor. The N=2 SUSY algebra implies the 
bound > 2|an + a^ml^ from below for their masses. Thus the states saturating 
this bound, known as Bogomolny-Prasad-Sommerfeld or BPS states, are described by the 
lattice spanned by the periods of dA. The singular locus of the fibration r(a. A), namely 
the points where the curve degenerates and a period aj or aD,j vanishes, corresponds then 
to vacua where one or several dyons become massless. 

For pure SU(2) Yang-Mills, the monodromy prescription at oo is restrictive enough to 
suggest the identification of the quantum moduli space of vacua with H/T{2) (H denotes 
the upper half space, and r(2) the subgroup of SL(2,Z) matrices congruent to 1 mod 2), 
assuming the minimal number two of singularities in the interior of the quantum moduli 
space. Since then, spectral curves have been proposed for a variety of gauge theories with 
matter, based on physical considerations such as decoupling, or analogies with singularity 
theory or soliton theory (see e.g. [40] and references therein). However, at the present 
time, we still do not have a complete correspondence between the group theoretic charac- 
terization of the gauge theory, consisting of the group G and the representation R for the 
hypermultiplets, and the fibration of spectral curves which characterizes its geometric and 
physical content. 

C. The Framework of the Theory of Solitons 

Nevertheless, an intriguing feature of most of the spectral curves for N=2 SUSY 
gauge theories known so far, is that they, together with the one-form dA, fit exactly in the 
framework of the foliation on Aig{n, m) with dA = QdE. In particular, 

• The spectral curves for SU(A^c) theories with Nf < 2Nc hypermultiplets in the 
fundamental representation of bare masses rrii, 1 < i < Nf, are given by the leaf (F, E, Q) 
with the following properties. 

- d-E has simple poles, at points P+, P_, Pi,with residues — A''^, N^. — Nf, and 1 (1 < 

i < Nf) respectively. Its periods around homology cycles are integer multiples of 27ri; 

* In this section, we adopt the normalization (6.9) for the periods of dA rather than 
(3.39), in keeping with the literature on Seiberg-Witten theory. Similarly, the present 
differs from the earlier r-function .Fwhitham of soliton theory (c.f. (4.35)) hy = 

27ri •^Whitham • 
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- (5 is a well-defined meromorphic function, with simple poles only at P+ and P_ ; 

- The other parameters of the leaf are fixed by the following normalization of the one- 
form dA = QdE 

Resp. (dA) = —nii 

Resp+(^dA) = -Nc2-^/^% Resp_(^dA) = {Nc - Nf){ )V(A^c-iv,) 

Resp^(dA) = (6.10) 

Here z = E~^l^'' or 2; = is as usual the holomorphic coordinate system 

near P+ or near P_ adapted to the Abelian integral E. 

These conditions imply that F is hyperelliptic, and admits an equation of the form 
(see [39]) 

= n - ^'^)' - A'^=-^^ + = ^(^)' - ^(^) (6.11) 
fe=i i=i 

Strictly speaking, the parameters afc of (6.11) agree with the classical vacua in (6.3) only 
when Nf. < Nf. For Nf > N^ there are 0(A) corrections, which can be absorbed in a 
reparametrization leaving the prepotential JF invariant [15]. Thus we may view the of 
(6.3) and (6.11) as identical. If we represent the Riemann surface (6.11) by a two-sheeted 
covering of the complex plane, then the meromorphic function Q on F in dA = QdE is just 
the coordinate in each sheet, while the Abelian integral E is given by £' = iog{y + A{Q)). 
The points P± correspond to the points at infinity, with the two possible choices of signs 
± for J/ = iV^^ — B. To choose a canonical homology basis, we let a;^, 1 < A; < Nc, 
be the branch points A{x^)'^ — B{x^), be a simple contour enclosing the slit from 
to for 2 < A; < A^c) cind B]^ be the curve going from to x'[ on each sheet. 
We can now give a preliminary and easy check that the curve (6.11) is consistent with 
the expected behavior of the theory at weak-coupling. Consider for simplicity the case 
of pure Yang-Mills, Nf = 0. Then as A — > 0, the discriminant of the curve behaves as 
-A^° Y\.j<ki'^j ~ ^fc)^7 s-iid the singularities are at the expected location. Furthermore, in 
this limit, dA ~ Q^^^dQ and the cycles A^ are just contours in the complex plane 

around a slit which shrinks to a single point a^. The residue formula gives at once 

afc = 2^ £ Q^^dg + 0(A^=) = afe + 0(A^O, (6.12) 

identifying ak as a classical order parameter. 

• The spectral curves for the other classical gauge groups with matter in the funda- 
mental representation are restrictions of the ones for SU(A^c) [4] [16]; 
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• The S'{](Nc) theory with matter in the adjoint representation is of particular interest. 
For massless matter, the theory has actually an N=4 supersymmetry, and is conformally 
invariant. As the hypermultiplet acquires mass, the N=4 SUSY is broken down to an 
N=2 SUSY. In [18], Donagi and Witten argued that the spectral curves for the theory 
are then given by Hitchin systems. Expressed in terms of elliptic Calogero-Moser systems, 
the curves they proposed are given precisely by the leaf (F, k, z) in Section III.D. Here the 
hypermultiplet mass has been scaled to 1, and the moduli r — 002/^1 of the torus is the 
microscopic gauge coupling. 

Although this suggests a deep relation between N=2 gauge theories and integrable 
models, such a relation is still not fully understood at the present time. Nevertheless, the 
parallelism between the two fields allows us to apply to the study of the prepotential of 
gauge theories the methods developed in the theory of solitons. Thus Theorem 16 implies 
readily [17] 

Theorem 20. The prepotential T for S\](Nc) gauge theories with Nf < 2Nc hypermul- 
tiplets of masses rrij in the fundamental representation, satisfies the following differential 
equation 

E «^ + E a;;^ - 2-^ = - ^ i^'^^ (.dA)Resp, (z-'dX) 
3=1 ^ 3=1 ^ 

+ Resp_(^dA)Resp_(z-MA)] (6.13) 

We observe that there is a slight abuse of language here, since in the case of the 
effective prepotential of gauge theories, is only fixed up to ak independent terms by 
(6.9). This is consistent with the fact that only derivatives of JF with respect to ak occur 
in the effective action. Thus the ak independent terms on the right hand side of (6.13) can 
be ignored by adjusting JF. The prepotential (4.35) (restricted to the leaf corresponding to 
the SU(A^c) theory) is one choice of J^. Another choice suggested by dimensional analysis 
(c.f. (6.8)) is the prepotential satisfying the homogeneity condition (^ + T>)J^ = 0. In 
this case, we recognize (6.13) as a renormalization group equation, with the beta function 
given by the right hand side of (6.13). Earlier versions of (6.13) appear in [24] [46] [47] [56] . 

To illustrate the power of Theorem 20, we shall show how it can generate explicit 
expressions for the contributions of instanton processes to any order. Thus we consider the 
regime where A is small and all the Ak-cjcles degenerate simultaneously. A fundamental 
observation is that in this regime, the quantum order parameters ak are perturbations of 
their classical counterparts which can be determined explicitly to any order. In fact, as 
noted in the arguments leading to (6.12), the Ak cycles are simple contours shrinking to a 
point in one sheet of the Riemann surface and residue formulae apply. The approximation 
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(6.12) can be improved to 



m=l ^ ' 



= iw^^ 

The evaluation of the dual periods aok is of course more difficult. We need to show that 
the prepotential JF, as defined by the S^-periods, reproduces the classical prepotential 
.7-"*-°^ in (6.8) (with hypermultiplets) and satisfies the non-renormalization theorem. This 
requires an analytic continuation in an aiixiliary parameter ^, as explained in [15]. How- 
ever, once this is established, the difficult instanton contributions can be derived from the 
renormalization group equation. Setting = :r(o) + ^(i) + :r(2) + • • ., we have, say up to 
2-instanton order and using Euler's homogeneity relation. 



(We note the overall factor Nf — 2Nc, which confirms the known conformal invariance of 
the theory with Nf = 2Nc- For the spectral curves of this theory, we refer to [4]). On 
the other hand, dropping all afc-independent terms, the right hand side of (6.13) is easily 
found 

i:''^s^+i:'»^*s--^^=«(^'-^^JE4. (6.16) 

3=1 ■> j=l J k=l 

We can now use (6.14) to reexpress the right hand side of (6.16) in terms of the quantum 
order parameters a^. The instanton contributions can then be read off after suitable 
rearrangements [15] [17] 

= 3^A''^"'-"'HE ^;l"'''tp'' +iE^^K)C5^K)) (6.17) 
-i^TTZ i^afc ai) 4^^^ 

Here the function Sk{x) is defined in analogy with (6.14) by Sk{x) — fTwr 

We turn next to a determination of the effective prepotential at strong coupling. In 
general, when a single cycle Ak or Bk degenerates, we expect the effective prepotential 
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to be expressible in terms of functions on the resulting surface of lower genus. A partic- 
ularly interesting case is the behavior of T near a point on the quantum moduli space 
of maximum degeneracy, where all B}, cycles degenerate simultaneously, and the spectral 
curve degenerates to two spheres connected by thin tubes. Physically, this means that a 
maximum number of mutually local dyons become simultaneously massless. As shown in 
[20], the points of maximum degeneracy occur at the curves = AQifff" — 4A^"'^=, where 
Ao(Q) is given by the A'"-th Chebyshev polynomial 

Ao = 2A^=Co(-^), Cq{z) = cos(iVarccos(^)) (6.18) 

A neighborhood of the maximum degeneracy point on the quantum moduli space is 
parametrized by polynomials P{Q) of degree N — 2, with the spectral curve = A{Q)'^ — 
4A2^'=, A{Q) = Ao{Q) + 2A^'=P(^). Since it is the Bk cycle which degenerates this time, 
it is more convenient to express the prepotential J-'{aDk) and the beta function in terms of 
the dual variables a^k, which can then be evaluated to an arbitrary order of accuracy by 
residue methods. The renormalization group equation remains the same under interchange 
of the dual variables ^ aok 

y]aDk4^-2J'= '^u{an) (6.19) 

where u is the coefficient of Q^=~^ in A{Q). Residue calculations show next aok is of first 
order in P and that 

oo , , 

aok = i{-fj^P{ck) + 2J «S ' ^fe = ^^^(]v")' = ^^^^Tj"^' ^ = 1' ■ ■ ■ ' (6.20) 

m=l ^ 

where the o)^^ are of order 0™'(a£)) and can be evaluated explicitly. To third order in 
o-Dki we find for u 

Solving the renormalization group equation, we obtain [14] 

Theorem 21. Near the 'point of maximum degeneracy on the quantum moduli space given 
by (6.19), the prepotential T is given by the following expression 

k=l k=l k=l 



Nc-1 2 

^ {ck-ci)^^ 



ki^l 
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up to third order in the order parameters a^k- Here is determined by log^ = |+logSfc. 

We should mention that there is by now an extensive hterature on Seiberg-Witten 
theories, and we refer to [40] for a description of other recent advances and for a more 
complete list of references. 
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